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Abstract 


This  report  details  the  data  and  program  structures  for  a  conver¬ 
sational  programming  system  which  translates  commands  describing  a 
Markovian  queueing  network  into  a  matrix  of  transition  intensities,  and 
which  provides  equilibrium  distributions  and  related  solutions  of  the 
network  according  to  requested  specifications. 
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Preface 


This  report  is  in  the  form  of  a  proposal  heavily  dosed  with 
tutoridl  matter.  It  was  origiruilly  written  in  the  latter  part  of  1967, 
before  work  had  begun  on  the  program  system  now  known  as  the  Queue 
Analyzer  System  (QAS).  A  working  version  of  QAS  [9]  was  completed 
ir  early  1969,  adhering  in  ail  major  respects  to  the  principles  laid 
focth  in  this  report,  except  that  the  solution  technique  of  Sections  6 
and  7  has  not  been  used  directly. 

While  the  system  proposed  here  is  described  in  terms  of  a 
specific  hardware  configuration,  this  work  is  not  limited  in  use  to 
any  one  such  configuration.  It  should  be  noted  that  for  a  program¬ 
ming  system  to  be  truly  efficient,  it  is  always  necessary  to  tailor 
its  structures  and  formats  in  some  way  to  the  hardware  used.  This 
has  been  done  here.  However  the  techniques  and  structures  are 
readily  adapted  to  other  cor  figurations. 

Because  of  the  unique  nature  of  the  system,  it  is  instructive 
tn  publish  this  report  in  essentially  its  original  form,  with  only  minor 
attempts  to  bring  it  to  correspondence  with  the  final  system.  .  hu?, 
it  is  more  a  treatment  of  programming  philosophy,  specialized  the¬ 
ory,  and  technique,  rather  than  a  complete  documentation  of  an 
existing  program.  Nevertheless-  all  essential  information  required 
to  understand  the  operation  of  QAS  is  here  contained,  as  well  as  some 
information  pointing  to  future  directions. 

xi 


1,  INTRODUCTION 


A  graphical,  problem-oriented  system  is  being  developed,  by 
meirs  of  which  solutions  to  simple  stochastic  networks  can  be  ob- 
ti’ped  rapidly  enough  to  facilitate  conversational  use.  The  problem 
descriptions  will  be  constructed  in  network  diagram  form  on  a  re¬ 
mote  DEC  339  graphic  console.  Simultaneously,  information  con¬ 
cerning  the  construction  will  be  sent  via  dataphone  to  a  time-shared 
IBM  360/67,  which  will  prepare  the  solutions.  It  is  the  purpose  of 
this  report  to  generally  specify  the  central  360/67  programs  and 
data  structures  which  accomplish  this  latter  feat.  These  programs 
will  be  .-eferred  to  here  as  the  Queue  Analyzer  System  (QAS). 

The  networks  drawn  will  be  restricted  by  the  pictorial  lang¬ 
uage  syntax  to  systems  which  can  be  modeled  by  a  continuous -time, 
f'r  'te  Markov  chain,  and  will  be  solved  by  numerical  solution  of  the 
Kolmogorov  equiiibr'um  equations.  The  advantage,  in  terms  of 
speed  and  precision,  of  this  technique  over  simulation  procedures 
’s  documented  in  [2]. 

Because  the  data  required  by  the  analysis  program  is  in  a  very 
different  form  from  that  describing  a  network  (which  is  in  terms  of 
blocks  and  connections),  a  translator  from  the  latter  to  the  former 
:s  required.  This  translator,  called  the  network  compiler,  is  the 


central,  and  most  difficult,  operation  carried  out  by  this  system. 
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'  In  addition,  there  are  four  other  major  operations  ^hich  the  system 

't  must  perform.  It  must  generate  the  network  description,  which  is 

I 

I  the  input  to  the  compiler,  from  information  contained  in  commands 

‘  supplied  by  the  console.  It  must  analyze  the  resulting  structure  for 

*  the  equilibrium  state  probabilities.  It  must  calculate  from  these  prob¬ 

abilities  the  specific  results  requested  by  the  console  user.  Finally, 
it  must  accept  new  definitions  of  symbols  used  in  network  generation. 

All  of  these  operations  are  to  be  performed  on  the  360/67.  The 
manipulation  of  pictures  (network  diagrams),  the  preparation  and 
updating  of  display  files,  and  the  recognition  of  commands  and  inter¬ 
rupts  from  the  user  are  all  solely  the  responsibility  of  the  SELMA 
system  in  the  remote  PDP9  which  is  part  of  the  333  console.  The 
SELMA  system  [ll]  must  also  keep  the  QAS  system  informed,  via 
the  dataphone  link,  of  any  operations  which  concern  it,  and  must  is¬ 
sue  commands  to  it  to  initiate  major  operatiorjs.  The  QAS  system 
will  carry  out  its  processes  under  the  control  of  the  remote  compu¬ 
ter,  treated  as  an  input  file  (probably  named  ^SOURCE*  in  MTS  [4] 
terminology).  It  will  send  its  responses  to  the  remote  computer  as 
an  output  file  (♦SINK*).  To  accomplish  this,  QAS  will  have  its  own 
version  of  the  current  "displayed  network"  stored  in  the  central  com¬ 
puter  memory. 

In  order  to  reduce  some  of  the  QAS  system-programming  prob¬ 
lems  to  manageable  proportions,  certain  limitations  of  ca^bility 
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have  been  accepted.  These  iimit-ations  are  chiefly  ones  which  limit 
the  meanings  which  can  be  assigned  to  network  symbols,  and  the  man¬ 
ner  in  which  the  symbols  can  be  related.  Our  minimal  objective  in 
this  system  has  been  to  provide  a  system  which  can  at  least  treat  net¬ 
works  consisting  wholly  of  queues,  exponential  servers,  infinite 
sources,  infinite  sinks,  random  branches,  merges,  and  priority 
branches.  While  many  other  symbols  can  also  be  treated  in  this  sys¬ 
tem  we  can  by  no  means  treat  dll  meaningful  symbols.  Nevertheless, 
those  which  can  be  treated  define  a  significant  class  of  models  having 
considerable  variety  and  power. 


2.  THE  BASIC  COMPILATION  AND  CAIXULATION  PHILOSOPHY 

A  continuous-tiTie.  finite  Markov  chain  is  a  stochastic  process 
which  takes  on  values  in  a  finite  set  of  points  callea  stales.  If  the  set 
of  states  is  mapped  by  a  one-one  function  h  onto  a  set  of  integers  |0. 1. 

- N  I .  and  if  the  equilibrium  probability  of  the  state  maopect  to  the 

integer  k  is  represented  by  tTj^.  then  the  vector  z  =  77^,.  . . . . 

^  s 

of  these  probabilities  is  a  solution  of  the  equation 

•.tU  =  0  (2.1) 

where  U  is  a  matrix  of  transition  intensities  descriptive  of  the  Maikov 
chain.  An  efficient  procedure  [l  j  for  solving  equation  (2. 1)  for  large 
N  has  been  in  use  for  several  years,  ana  derives  much  of  its  efficiency 
from  the  form  in  which  U  is  stored.  An  improvea  storage  structure  for 
this  purpose  calleo  the  matrix  outline  structure  [Sj  has  been  devisea. 
and  proposed  for  this  project.  An  adaptation  of  the  earlier  procedure 
to  this  new  structure  will  be  required,  and  will  be  called  RQA-?.  Its 
purpose  is  simply  to  calculate  a  when  U  is  given. 

In  contradistinction  to  equation  (2. 1).  the  description  of  the  Mar¬ 
kovian  system  provided  by  the  console  is  in  the  form  of  a  network.  We 
may  say.  by  way  of  definition,  that  a  network  N  consists  of  a  set  .. 
of  elements  and  a  set  ^  of  connections. 

N  -  :  .  -  (2.2) 


I 


I 

{ 


j  • 
|- 
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where  the  terms  element  and  connection  are  yet  to  be  defined  but  have 
the  intuitive  meaning  implied  by  "blocks  having  a  distinctive  function" 
and  "lines  joining  the  blocks  in  particular  ways. "  Our  objective  is  the 
construction  of  tlie  matrix  outline  structure  for  the  matrix  U  from 
knowledge  of  the  structure  that  describes  the  network  N.  It  is  this 
process  which  we  term  compilation  of  the  network. 

It  should  be  noted  that  the  state  of  the  Markov  chain  will  be  re¬ 
lated  to  the  joint  conditions  of  all  of  the  elements  in  the  network,  so 
that  a  part  of  the  compilation  will  be  to  form  the  state  mapping  from 
a  set  of  multi- dimensional  state  vectors  to  the  set  of  consecutive  inte¬ 
gers  {0, 1, 2, . . . ,  N  }.  It  is  only  after  this  mapping  is  specified  that 
s 

the  matrix  U  defines  a  Markov  chain  unambiguously.  (The  process  of 
deriving  a  suitable  mapping  is  not  a  trivial  one,  since  the  boundaries 
of  the  set  of  possible  states  will  generally  be  irregular,  and  the  map¬ 
ping  must  be  one- one  on  consecutive  integers  if  computer  memory  is  to 
be  used  efficiently. ) 

The  calculation  process  witli  which  we  are  concerned  involves 
eight  distinct  phases. 

(1)  Generation  of  the  Network 

(2)  Compilation  of  the  Network 

(3)  Definition  of  New  Elements 

(4)  Calculation  of  Equilibrium  Probabilities 

(5)  Specification  of  Results 
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(C)  Calculation  of  Results 

(7)  Documentation 

(8)  System  Control  and  Supei'vision 

The  user  at  the  console  creates  a  network  through  a  long  series 
of  actions,  such  as  creation  of  elements  or  connections,  evaluating 
parameters,  changing  connections,  etc.  The  sum  total  of  all  these 
commands  defines  the  network.  The  SELMA  system  (or  possibly  a 
teletypewriter)  conveys  this  list  of  commands  to  QAS,  which  must  gen¬ 
erate  the  appropriate  descriptive  structure,  step  by  step.  This  process 
is  the  process  of  network  generation  (phase  1,  above).  In  definition  of 
new  elements,  a  series  of  commands  is  received  which  represents  the 
meaning  of  a  new  symbol.  This  nieani^'g  must  be  assimilated  during  a 
definition  phase  (phase  3),  and  preserved  in  a  form  which  can  be  used 
to  identify  future  occurrences  of  the  new  symbol.  Similarly,  when  de¬ 
sired  results  are  being  specified  '.phase  5\  the  specifications  are  re¬ 
ceived  as  a  stream  of  commands  whose  meaning  must  be  collected  in 
a  specific  form.  This  form  is  u*jed.  by  the  calculator  program  (phase 
6)  to  prepare  a  matrix  which  can  jjroduce  the  desired  result  when  ap¬ 
plied  to  the  vector  of  equilibrium  state  probabilities.  The  documenta» 
tion  phase  (phase  7)  is  one  in  which  names  are  associated  with  struc¬ 
tures  to  be  saved,  and  structures  in  the  system  can  be  substituted  by  a 


named  file. 
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Each  of  the  phases  is  a  distinctly  different  operation,  requiring 
its  own  set  of  programs  and  data  structures.  Each  Is  likely  to  proceed 
for  a  significant  interval  of  time,  and  to  continue  at  the  discretion  of 
the  SELMA-generated  commands  which  are,  in  turn,  responsive  to  the 
user's  char;ges  of  pace  or  purpose.  Any  phase  can  be  terminated  be¬ 
fore  completion  and  resumed  later.  It  is  the  purpose  of  the  "System 
Control  and  Supervisor"  programs  (phase  8)  to  recognize  changes  of 
phase  from  the  information  received  from  SELMA,  and  make  appropri¬ 
ate  changes  in  the  current  files  contained  in  (virtual)  memory.  (Since 
the  entire  system  is  large,  and  the  user's  sessions  are  likely  to  be 
long,  it  would  be  uneconomical  to  keep  all  files  resident  throughout 
the  session.)  The  segmented  nature  of  the  360/67  addressing  scheme 
seems  particularly  well  suited  to  this  kind  of  overlaying  of  programs. 

In  the  succeeding  sections  of  this  report,  the  structu^'es  of  the 
special  files  required  will  be  described,  along  with  the  programs  Ln  the 
various  phases  which  use  them.  In  the  remainder  of  this  serttor.  some 
further  comments  about  the  Supervisor  and  general  file  crgarization 
will  be  made.  The  Supervisor  program  is  tentatively  diagrammed  in 
Figure  2. 1. 

The  Supervisor  will  be  assumed  to  reside  in  virtual  memory 
throughout  the  session.  One  of  its  chief  component  parts  is  a  program 
which  regularly  reads  the  commands  and  data  arriving  from  SELMA. 
These  are  buffered  (queued),  and  treated  in  the  order  in  which  they 
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arrived.  Thus,  QAS  operates  asynchronously  with  SEL>4A,  lagging  be¬ 
hind  it  more  or  less  depending  upon  360  scheduling,  rat^  of  SELMA  com¬ 
mand  generation,  and  its  own  speed.  For  this  reason,  parring  blunders 
of  either  system,  care  must  be  exercised  that  no  comtn^ad  received 
from  SELMA  and  queued  by  the  Supervisor  should  be  able  to  produce  a 
"fatal"  type  error  in  QAS.  QAS  must  not  grind  to  a  halt  When  SELMA 
commits  a  syntax  error,  for  example,  because  by  the  time  the  error 
flag  is  conveyed  to  SELMA,  SELMA  may  have  gone  far  DoyOnd  the  point 
of  error  and  be  unable  to  retrace.  The  result  would  be  a  dilemma  re¬ 
solvable  only  by  starting  again  from  a  clean  slate — an  intolerable  burden 
for  the  user.  On  the  other  hand,  after  certain  commands  3ELMA  might 
wait  for  an  affirmative  response  before  issuing  more.  One  such  case 
would  be  after  the  command  for  a  network  compilation.  Ka  other  oper¬ 
ations  should  proceed  until  it  is  known  the  network  was  confipilable  and 
did  not  require  changes. 

The  Supervisor  has  responsibility  for  a  table  of  af'sa  names,  and 
for  calling  the  routines  for  each  particular  phase  via  a  LlhJlC-type  call 
to  MTS.  Every  area  should  have  its  own  table  of  special  pointers  to 
items  within  itself,  the  structure  of  the  table  being  kno'Wu  to  every  pro¬ 
gram  which  uses  it.  In  that  way,  the  Supervisor  does  pot  peed  to  keep 
tables  other  than  the  area  tables. 


$ 


3.  THE  GENERATION  PROCESS  AND  NETWORK  SYNTAX 

3.1  Elements  and  Connections 

An  element  is  defined  by  its  type,  its  parameter  value.s,  and  its 
ports.  The  type  r  of  element  e  describes  the  nature  of  the  element. 
Examples  of  types  of  elements  are  "queue."  "server."  "exit."  and 
"source.  "  The  parameter  values  modify  the  type  for  the  particular 
element.  Thus,  in  the  case  of  the  queue  there  is  one  parameter,  its 
capacity.  For  the  server,  by  convention,  it  is  the  reciprocal  of  the 
mean  holding  time.  The  exit  and  source  have  no  parameters.  Other 
elements  may  exist  which  have  more  than  one  parameter.  Tnus.  for 
each  element  the  parameter  values  can  be  expressed  in  general  as  a 
vector,  or  list.  The  number  of  parameters  of  an  element  is  a  prop¬ 
erty  which  is  characteristic  of  the  type  of  Uie  element.  A  port  is  an 
ab.stract  object  which  represents  a  "place  for  joining"  elements  to 
other  elements  via  connections.  The  queue  and  server  each  have  two 
ports.  The  exit  and  source  have  one.  The  number  of  ports  are  thus 
clearly  also  a  characteristic  of  the  type  of  the  element. 

Let  the  set  of  elements  v'  of  the  network  be  an  indexed  set 
over  some  set  I : 

O  =  {e.:  ie  I|,  (3.1) 

Let  r.  be  the  type  of  the  element  e-  letp.  be  the  sequence  of  its 
parameter  values,  and  let  be  a  set  of  ports. 


11 


12 


Pj  --  Ipjj;  i  6  Jj!.  (3.2) 

Then  the  element  e.  can  be  al<iebraically  describcfi  in  the  form 

e.  -  -  .p..  P.  '  .  for  all  i  f  I  (3.  3) 

For  certain  useful  elements,  although  not  for  tlie  four  basic  primitives, 
the  number  of  parameters  and  the  number  of  ports  in  a  particular  ele¬ 
ment  is  a  function  of  one  or  more  of  the  parameters.  Since  tne  nature 
of  Pj  and  P.  is  consequently  dependent  upon  tne  value  of  this  parameter 
(or  parameters),  this  parameter  (or  jxirameters)  must  be  supplieo  be¬ 
fore  equation  (3.3)  has  meaning.  Such  jxirameter(s)  will  be  called  gen- 
era t  ion  parameters  because  of  tneir  special  role  in  the  network  genera¬ 
tion  phase.  They  will  always  be  placed  first  in  the  p;irameter  list. 

Quite  similarly  to  elements,  connections  are  defined  by  their 
type,  their  parameter  values,  and  the  ports  they  connect.  Types  t)f 
connections  in  our  basic  list  are  '’simple.”  "overflow  bran<“!i”  fsiime- 
times  callea  "priority  branch"),  "join.”  ana  "random  branch."  The 
remarks  concerning  jxiranieter  values  of  elements  apply  equally  well 
to  those  of  connect  ions.  The  p->rts  connected  by  connections  are  the 
ixjrls  of  the  elements,  and  each  port  may  be  involved  in  only  one  con¬ 
nection.  Thus,  if  the  connection  set  is  represented  by  the  inae.xed 
set 

icj.:  k.  Kl.  (3.4) 

then  a  connection  Cj^  is  reprc.m'ntcd  as 
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(3.5) 


where  and  Pj^  are  the  type  and  the  parameter  list,  respecth'eiy,  and 
is  a  subset  of  tne  ports  of  the  elements 

P.  C  -  Pj  (3.6) 

^  ~  i  e  I  ' 

Every  port  is  in  exactly  one  connection: 


Pk  Pk  =  <«>  kj.k^eK  (3.7) 

P,.  =  P..  (3.8) 

k  €  K  ^  i  £  I  ' 

(In  other  woras  |Pk*  ^  Kj  ana  jP. ;  i  c  l]  are  each  partitionings  of 
the  set  of  ail  ports. ) 

The  random  branch  is  an  example  of  a  connection  having  a  gen¬ 
eration  parameter.  This  parameter  is  the  number  of  "branches” 
among  which  tasks  can  be  switched.  The  number  of  ports  connected  by 
this  connection  is  one  more  than  the  number  of  branches,  as  is  the 
number  of  parameters. 

The  elements  and  connections  as  described  here  collectively  aes- 
cribe  a  netv/ork.  However,  it  is  not  necessary  that  these  elements 
correspond  precisely  with  the  symbols  used  at  the  console.  This  no¬ 
tation  represents  a  language  for  describing  networks  wnich  is  distinct 
from  that  usea  by  SELMA,  ana  it  is  the  PDP-9’s  responsibility  t\)  trans¬ 
late  between  these  languages  when  issuing  generation  commands  to 
QAS.  Nevertheless,  if  we  associate  particular  graphs  with  each  element 
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type  in  QAS.  then  n  netwox*k  N  can  be  drawn  like  that  in  Figure  3. 1. 


Figure  3. 1 . 

A  Diagram 

Such  a  diagram  is  useful  in  visualizing  networks,  and  may  be  so  used 
in  this  report.  A  summary  of  properties  of  the  standard  elements 
and  connection  types  is  given  in  Table  3. 1.  along  with  a  set  of  graph¬ 
ics  which  may  be  used  for  them. 

3. 2  Syntax  of  Networks  and  Diagrams 

Equations  2. 1  and  3. 1  through  3.  8  define  the  syntax  of  a  net¬ 
work  as  required  by  QAS.  QAS  must  enforce  rules  upon  the  specifi¬ 
cation  of  networks  to  it  which  will  guarantee  that  the  netwc)rk  satis¬ 
fies  this  form.  In  addition.  SELMA  may  apply  further  rules  for  the 
form  of  its  diagrams.  These  are  primarily  to  protect  the  user  from 


Table  3. 1  Connection  Characteristics  and  Symbols  for  Standard  Objects 
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constructions  which  may  be  compilable,  but  not  meanifigfui  in  semantic 
terms. 

The  chief  syntactic  rule  of  QAS  is  that  ports  of  elements  may  only 
be  joined  through  connections.  A  second  constraint  is  the  one,  mentioned 
in  the  description  of  connections,  that  requires  ports  to  be  associated 
with  exactly  one  element  and  exactly  one  connection.  A  third  constraint 
is  that  parameters  of  an  object  (element  or  connection)  must  be  numer¬ 
ically  valued  (not  variables)  at  the  time  of  compilation,  and  generation 
parameters  must  be  numerically  valued  at  the  time  of  the  fir-st  mention 
of  the  object  (i.  e. ,  at  the  time  of  creation).  Finally,  elements  and  con¬ 
nections  are  restricted  to  instances  of  types  whose  meaning  has  been 
previously  defined  using  the  notation  to  be  described  in  sections  4  and 
5  of  this  report.  These  restrictions  merely  reflect  the  requirements 
necessary  to  fit  the  data  structures  used  by  QAS. 

The  first  constraint,  part  of  the  second  constraint  (the  require¬ 
ment  that  exactly  one  element  be  associated  with  a  port),  and  the  last 
constraint  will  all  be  automatically  enforced  because,  using  the  com¬ 
mand  set  to  be  described,  SELMA  will  be  unable  to  violate  them.  How¬ 
ever,  it  is  possible  to  attempt  multiple  connections  at  a  port,  or  to 
leave  a  port  unconnected,  or  to  supply  the  wrong  number  of  genera¬ 
tion  parameters,  or  to  fail  to  provide  values  for  other  parameters. 

Either  SELMA  will  be  expected  to  guarantee  that  the  proper  constraints 
are  met,  so  that  QAS  may  treat  their  violation  as  fatal,  or  QAS  will  be 
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empowered  to  test  for  the.R  at  the  time  tiieir  violation  can  cause  diffi¬ 
culty.  and  demand  immediate  correction  from  SELMA.  The  latter 
course  means  that  there  must  be  a  procedure  for  the  QAS  supervisor 
(the  only  program  residing  in  memory  at  all  times)  to  interrupt 
SELMA  demanding  a  specific  type  of  response  without  losing  commands 
from  the  input  stream.  In  any  case,  there  must  be  provision  in  QAS 
to  test  these  properties  at  appropriate  times,  and  to  return  diagnostics. 

SELMA  will  require  that  ports  be  characterized  as  either  input 
or  output,  and  that  each  connection  type  may  connect  inputs  and  outputs 
only  in  certain  ways.  QAS  does  not  have  this  requirement,  and  will 
not  test  for  it.  Indeed,  it  will  not  even  know  which  ports  are  inputs 
and  which  are  outputs. 


auiw  ui 
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<  n  , . . .  > 
^1 


P2 

P. 


Po  = 


<  n  , . . .  > 
^2 

{ P^’  Pg*  •  •  •  } 

{pg> . . .  }. 


Figure  3.2  describes,  schematically,  a  data  structure  for  this  network. 
The  meaning  of  the  various  diagram  symbols  is  called  out  in  Figure  3. 3. 
We  have  represented  sets  by  a  one-word  head  of  a  ring,  the  ring  linking 
its  elements.  The  values  of  n-tuples  have  been  represented  by  an  n- 
word  contiguous  block,  while  elements  belonging  to  a  predeterminable 
number  of  sets  have  been  represented  by  a  like  number  of  contiguous, 
one-word  ring  elements  followed  contiguously  by  the  value  of  the  ele¬ 
ment. 


Because  every  set  in  this  structure  contains  elements  having  iden¬ 
tical  form  and  because  the  programs  operating  on  these  structures  will 
always  know  what  kind  of  element  is  being  operated  upon,  the  use  of  con¬ 
tiguous  blocks  in  this  manner  is  feasible,  and  will  result  in  an  efficient 
and  compact  structure.  If  objects  belonging  to  a  variable  number  of 
sets  were  present,  a  form  of  association  structure  such  as  that  used 
by  Newman  [5]  could  be  easily  incorporated  using  "nubs”  and  nested 
rings  in  addition  to  the  above.  So  far,  that  has  been  unnecessary. 
Although  the  structure  of  the  network  is  to  be  manipulated  in  the  course 
of  executing  these  programs,  the  sizes  of  contiguous  blocks  will  not 
change  during  the  execution.  However,  blocks  will  be  created  and 


A  fragment  of  a  completed  network  structure. 
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Ring  Element  (Defining  Set) 
Ring  Element  (Definition) 
Value 

Dummy  Value 
Pointer 


-<nj 

-HZZl- 

IZID 

CO 


Figure  3.3 

Network  structure  diagram  notation. 
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destroyed,  and  this  structural  technique  depends  heavily  upon  the  avail- 
ability  of  an  efficient  free-core  manager.  The  manager  proposed  in  a 
separate  memo  [6]  (Appendix  A)  appears  to  satisfy  the  requirements. 

3.4  The  ’’Generate"  Program 

The  generate  phase  program  consists  of  a  collection  of  routines 
corresponding  to  the  generation-commands  receivable  from  SELMA. 
Upon  entry  from  the  Supervisor,  the  generate  program  examines  the 
command  and  issues  a  call  to  the  appropriate  generation  routine.  The 
network  description  is  rec;eived  by  QAS  in  the  form  of  a  stream  of 
these  commands.  The  functions  of  the  "generation  routines"  carrying 
out  the  commands  are  summarized  in  Table  3. 2.  More  generation 
routines  may  be  added  as  use  of  this  system  oecomes  more  sophisti¬ 
cated. 

These  routines  operate  entirely  on  an  area  called  the  network 
area  of  meuiory  which  contains  the  network  structure,  the  type  struc¬ 
ture  (to  be  described  in  section  4),  and  a  symbol  table.  The  symbol 
table  provides  the  means  for  converting  element  names,  coimection 
names,  and  other  names  used  SELMA  to  corresponding  addresses 
in  the  network  area.  The  table  grows  as  commands  are  received  de¬ 
fining  new  objects.  Objects  also  can  be  removed  from  the  symbol 
table.  The  simplest  form  for  it  is  obtained  if  SELMA  assigns  consec¬ 
utive  integers  as  names  for  objects,  whereby  the  table  consists  solely 
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of  a  vector  of  pointers.  The  location  of  this  vector  within  the  network 
area  should  be  permanent.  The  symbol  table  should  also  keep  a  pointer 
to  the  network  block  (N)  as,  say,  its  first  symbol.  Thus,  the  symbol 
table  will  keep  the  locations  of  all  master  objects  as  well  as  those 
which  must  be  referred  to  by  SELMA.  This  will  permit  master  blocks 
to  be  moved  at  will  for  such  things  as  documentation  functions  which  use 
the  network  area,  by  merely  changing  a  pointer  in  the  network  area's 
symbol  table. 

To  facilitate  discovery  of  unconnected  ports,  which  would  prevent 
execution  of  a  compile  command,  a  special  master  object  R  iiS  in  tte 
symbol  table  representing  the  set  of  unconnected  ports.  When  elements 
are  created,  their  ports  are  automatically  joined  to  this  set.  "Connect" 
commands  unjoin  ports  from  this  set,  and  join  them  to  a  newly  created 
connection.  The  symbol  table  and  unconnected  port  set  block  are  illus¬ 
trated  in  Figure  3. 4. 

The  '  create  element"  and  "alter  generation  jarameter"  routines 
will  be  described  in  detail  in  section  4,  which  treats  the  information 
needed  to  define  types.  The  other  generation  routines  in  Table  3. 2 
should  be  self-explanatory. 

3  5  Service  Programs  Required 

Four  sets  of  service  programs  will  be  needed  by  the  generation 
routines.  They  should  be  FORTRAN-cailable  so  that  the  generation 
subroutines  can  be  written  in  FORTRAN.  They  will  also  be  found 
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useful  to  other  phase  programs,  such  as  the  compiler.  Their  func¬ 
tions  and  calls  are  summarized  in  Table  3. 3. 

The  list  of  programs  shown  is  probably  not  complete,  and  may  be 
added  to  once  programming  is  begun.  In  particular,  the  ring  operators 
represent  a  subset  of  the  commands  used  in  Roberts'  CORAL  system  [7], 
and  some  of  the  other  CORAL  commands  may  be  also  desirable. 
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Table  3.  3  Service  Programs  for  Generation — Summary 
(continued) 


4.  PRIMITIVE  ELEMENTS  AND  THEIR  MEANING 


The  elements  in  the  network  structure  have  been  described  in  terms  ; 

! 

of  their  types.  Examples  cited  were  the  types  "queue, "  "server,  "  "exit, "  | 

I 

5 

and  "source. "  Each  of  these  are  called  primitive  types  because  they  I 

{ 

i 

cannot  be  described  in  terms  of  any  other  element  types.  However,  a  \ 

description  oi  theic  meaning  must  be  given  if  the  compiler  is  to  have  suf¬ 
ficient  information  on  which  to  operate.  In  this  section,  a  means  for  for-  f 

mal  description  of  the  meaning  of  a  type  will  be  given,  and  a  data  structure 
for  conveying  that  meaning  to  the  translator  will  be  provided. 

4. 1  The  State  Set  , 

Every  element  has  a  set  of  states  associated  with  it.  These  states  i 

represent  the  various  possible  identifiable  conditions  which  the  element 
can  assume.  In  QAS,  the  state  set  of  any  element  will  be  a  finite  set  of 
non-  negative  integer  n-tuples,  where  n  is  characteristic  of  the  element. 

In  the  case  of  a  queue  whose  "maximum  length"  parameter  is  N,  this 

set  is  the  one-dimensional  set  of  integers,  {0,1 . nJ,  representing 

"the  number  of  waiting  tasks. "  In  the  case  of  the  server  it  is  the  one¬ 
dimensional  set  of  integers  {0, 1.  ?}  representing  the  "idle, "  "busy, " 
and  "holding"  caiditions.  respectively.  Finally,  in  the  case  of  both 
the  exit  and  source  it  is  the  set  {o}  a  condition  representing  "readi¬ 
ness  to  receive  (or  sendl  a  task. "  Although  all  four  of  the  basic  ele¬ 
ment  types  have  one- dimensional  state  sets,  higher  dimensions  are 
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possible  for  yet  unnamed  elements. 

The  things  which  most  characterize  elements  are  the  stimuli  which 
can  produce  changes  in  state,  and  the  responsive  changes.  These  are 
characterized  by  the  so-called  events  of  the  element-type.  Some  of  the 
stimuli  are  self-generated,  as,  for  example,  the  "service  completions" 
of  a  server.  Others  are  external  stimuli,  as  the  occurrence  of  "inputs 
of  tasks"  at  a  port  of  an  element.  The  stimuli  and  responses  of  the  for¬ 
mer  will  be  described  by  autoevents,  while  those  of  the  latter  will  be 
described  by  exoevents. 

4.2  Autoevents 

By  definition,  an  autoevent  of  an  element  e.  having  a  state  set 
Sj  is  a  triple 

where : 

b-  is  a  subset  of  S. 

Jr  1 

g^  is  a  constant  n-tupie,  such  that  for  each  x  £  b^, ,  x  +  g^ 
is  in  S. 

is  a  positive  real  number 
I  is  an  index  in  some  index  set  L.. 

Here  b^  is  called  the  autocondition  set  of  the  event  ,  and  represents 
a  set  of  states  for  which  the  event  is  oossible;  the  constant  u.  is  called 
the  rate  of  the  event,  and  represents  the  probability  intensity  of  its 
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occurrence;  the  constant  n-tuple  is  called  the  increment  of  the  event, 
and  represents  the  change  in  state  which  results  from  the  occurrence  of 
the  event. 

To  illustrate,  consider  the  "service  completion"  event  for  a 
server  having  mean  service  rate  (parameter)  y.  The  event  can  occur 
only  i:.  .ne  busy  state,  state  1,  and  w’hen  in  that  state  it  occurs  with 
orobability  intensity  y.  The  event  causes  a  change  from  state  1  to  the 
holding  state,  state  2,  so  that  the  change  in  state  is  unity.  Thus, 

t  .<{1},  1,  y>.  (4.2) 

It  is  possible  for  an  element  to  have  more  than  one  autoevent. 
Consequently,  for  every  element  there  is  a  set  of  autoevents,  called 
the  autoevent  set 

H.  -  (fjp  f  €  L,}  (4.3) 

of  the  element  e.  where  the  index  sets  L.  are  disjoint  over  all  i  in  I. 
Frequently,  when  the  format  of  equation  4. 1  is  insufficiently  general 
to  describe  a  particular  stimulus  and  its  response,  the  desired  des¬ 
cription  can  be  obtained  b>'  splitting  the  description  up  into  several 
autoevents. 

4. 3  Exoevents 

By  definition  an  exoevent  at  a  port  p  of  an  element  e.  hav¬ 
ing  a  state  set  S.  is  also  a  triple 


where: 


b_  is  a  subset  of  S. 
tn  1 

g  is  a  constant  n-tuple  such  that  for  each  x  €  b  , 
xtg^isinSj 

is  a  probability 

m  is  an  index  in  some  index  set. 

Here  b  is  the  endocondition  set  of  the  event,  and  represents  a  set 
m  - - 

of  states  for  which  the  event  is  possible  (e.  g.,  if  the  port  p  is  an  input 
poll,  b  is  the  set  of  states  for  which  an  input  of  a  task  can  be  accepted); 
the  probability  of  the  event  is  the  probability,  given  that  the  stimu¬ 
lus  (e.  g. ,  an  input  of  a  task)  has  occurred  and  that  the  element  is  Ln  a 
state  of  b^  that  the  change  in  state  produced  is  equal  to  g^,  the  incre¬ 
ment  of  the  event. 

An  example  of  an  exoevent  is  the  "input  of  a  task"  event  (call  it 
^  )  at  the  input  port  p  of  a  queue  whose  maximum  length  is  N.  An 

input  can  be  allowed  only  when  the  state  is  less  iMn  N,  and  it  results, 
with  certainty,  in  an  increase  of  the  state  by  unity.  Thus 

C  =  <  1,  1  >,  (4.5) 

for  N  >  1. 

For  every  port  p  of  an  element  e.,  there  is  an  exoevent  set  Z  ^(p). 
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representing  all  the  possible  responses  to  the  external  stimulus  at  the 
port  p.  The  function  Z  .  will  be  called  the  exoevent  function  of  the  ele¬ 
ment  e..  The  exoevent  set  is  a  set  of  exoevents  whose  probabilities 
add  up  to  zero  or  one  for  each  state  of  the  element.  Let 

Zi(p)  -  me  M.(p)},  p=  P^,  i  e  I,  (4.8) 

and 


’m 


b  ,  g  ,  TT  >  , 
m  ^m  m 


(4.7) 


where  the  b  ,  g  ,  ff  have  the  forms  appropriate  to  exoevents  of  an 
element,  and  where  the  M.(p)  are  disjoint  for  all  p  aud  i.  Then 


^  If  -  0  or  1  (4. 8) 

m: 

m 

for  each  x  e  S.  where  S.  is  the  state  set  of  the  element.  Those  states 
for  which  this  sura  is  zero  are  states  for  which  the  external  stimulus 
cannot  occur,  such  as,  for  example,  the  states  where  an  input  task 
cannot  be  accepted. 


4,  4  Type  Definition 

The  meaning  of  the  ’’type”  of  a  particular  element  e^  is  con¬ 
sequently  conveyed  by  the  triple  <  S^,  H.,  Z  .  >  ,  so  that  we  can  write 


r.  -■  Zj> 


(4.9) 


for  all  i  c  I.  The  above  rules  for  determining  the  state  set,  autoevent 
set,  and  exoevent  function  have  been  applied  for  instances  of  each  of 
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the  four  basic  element-types,  and  these  properties  have  been  sum¬ 
marised  in  Table  4.1.  In  each  of  the  first  two  cases,  the  port  p^  is 
the  input  port,  and  port  P2  is  the  output  port.  It  is  interesting  to  note 
that  the  exit  and  source  elements  are  identical  in  these  terms.  Con¬ 
sequently,  they  will  both  be  treated  as  instances  of  a  single  element- 
type,  the  source-exit.  Other  element-types  are ,  of  course,  possible 
provided  their  stimuli  and  responses  can  be  completely  described  by 
events  of  the  form  described. 

The  information  conveyed  in  Table  4. 1,  and  extensions  of  it  if 
more  element-tyi'es  are  defined,  must  be  made  available  to  the  com¬ 
piler.  In  particular,  there  must  be  a  data  structure  which  describes 
the  type  r.  of  an  element  e.  in  literal  form,  with  all  dependencies 
upon  parameters  evaluated  and  readily  available.  In  addition,  the 
compiler  in  the  course  of  its  v/ork  will  be  creating  more  elements 
which  are  not  instances  of  these  basic  types  but  are  eo'jally  described 
by  triples  <  S;,  H,,  >  .  The  structure  for  providing  this  informa¬ 

tion  to  (he  compiler  is  added,  at  the  appropriate  times,  to  the  network 
structure  by  use  of  the  r-line  in  the  element  blocks.  This  structure 
is  illustrated  schematically  in  Figure  4. 1  for  an  element  e^  such  that 


Pj  --  iPyPzl 

(4.10) 

(4.11) 

ZjCPi)  = 

(4.12) 

Z,(P2)  "  {Cj}- 

(4.13) 

Fragment  of 
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The  sets  of  states  Sj,  bj,  b^,  . . . ,  bg,  and  the  increments  g^,  . . . ,  gg 
are  not  shown;  their  form  will  be  discussed  in  section  4. 6.  The  block 
-j  'S  referred  to  as  a  type  block. 

It  should  be  noted,  however,  that  at  the  time  the  element  block 
was  created  (during  generation  phase)  the  type  line  rin  the  element  block 
contained  a  name  of  the  type,  rather  than  a  pointer  to  a  structure  like  that 
shown  in  Figure  4. 1,  It  is  not  possible  during  element  creation  to  form 
this  type  structure  because  the  parameter  values  are  not  yet  known.  (It 
should  be  noted  in  Table  4. 1  that  variables  N  and  y  were  necessary  to 
represent  the  state  set  and  the  events. )  It  is  also  not  completely  de¬ 
sirable  to  supply  this  structure  immediately,  even  if  it  were  possible, 
because  it  would  considerably  increase  the  amount  of  storage  required 
during  generation.  This  information  is  needed  only  during  compilation, 
or  if  compilation  has  been  partially  performed. 

Ir  any  case,  the  t  -line  in  the  element  block  has  two  possible  mean¬ 
ings,  and  during  compilation  either  meaning  may  appear  in  the  line. 

Since  a  pointer  requires  only  24  bits,  and  the  name  (a  type  number, 
actually)  can  be  contained  in  an  8-b't  byte  or  less,  it  is  not  difficult  to 
distinguish  between  the  two  cases.  By  reserving  the  name  ”00’’  (hex¬ 
adecimal)  to  literal  valued  types  (i.  e.,  those  for  which  a  structure  is 
present),  and  using  the  format  shown  in  Figure  4  2a  simple  routine 
can  determine  whether  name  or  pointer  is  present.  If  the  pointer  to  the 
type  structure  is  present,  we  say  the  type  is  literal  valued.  Otherwise, 


the  type  is  said  to  be  namc<l 


NAME 


POINTER 


0  78 
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Figure  4.  2  Format  of  r-line  of  element. 

4.  5  Type  Structure  Generation 

Wb«;.  ui€  compiler  creates  an  element  itself,  it  will  always  have 
a  literal  type.  When  the  element  is  created  by  the  "create  element" 
routine  during  generation  phase,  it  will  have  a  named  type.  Tht  con¬ 
version  of  the  named  type  to  a  literal  type  will  be  carried  out  only  when 
the  type  structure  of  the  particular  element  is  needed  by  the  compiler. 
The  compiler  requests  the  literal  type  of  an  element  through  a  routine 
called  the  type  finder.  The  element  is  specified  and  if  tne  element  is 
already  of  literal  type,  the  pointer  to  the  type  block  is  in  mediately  re- 
tur.'.ed.  If  the  element  has  a  named  type,  the  finaer  calls  a  particular 
type-evaluation  routine  which  create^  the  type  structure,  using  the 
parameter  values  of  the  element.  There  is  one  type- evaluation  routine 
for  each  type  name.  These  routines  are  written  wnenever  a  new  type 
name  is  added  to  the  repertoire  ofQAS. 

The  operation  of  the  type  structure  programs  is  summarized  in 
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Table  4.  2.  These  programs  are  used  during  compile  phase,  and  re¬ 
quire  the  presence  of  the  working  network  structure  only. 

4.  6  Representation  of  Sets  of  States  and  Increments 

The  £  et  S  for  each  element  is  a  set  of  n-tuples.  The  autocondition 
sets  and  endoccndition-sets  of  the  everts  of  the  element  are  subsets  of 
this  set.  For  reasons  of  efficiency  in  the  operation  of  the  compiler,  it 
must  be  possible  to  rap’dly  form  aniors.  intersections,  and  differences 
of  these  sets  For  efficient  use  of  storage,  these  sets  cannot  be  repre¬ 
sented  by  rings  of  members .  their  number  is  too  great. 

For  these  reasons,  a  spec’ll  data  structure  must  be  used  for 
represertirg  sets  of  r -tuples.  The  basic  unit  in  this  structure  will  be  a 
structure  representing  rectangular  svis  of  n-tuples,  sets  which  are  Car 
tesian  products  of  sets  of  consecutive  ’rtegers.  Any  state  set,  or  sub¬ 
set,  will  be  represerted  by  a  ur’cr  of  rectangular  sets.  This  structure 
is  illustrated  in  F'gure  4.  3.  The  basic  lire  describing  the  set  is  a  sin¬ 
gle  word  whose  format  is  shown  Figure  4.  4.  This  word  contains  the 
dimension,  n,  of  the  n-^uples  and  is  a  ring  head.  The  upper  byte  of 
the  rmg  elements  pointed  tn  is  zero,  ard,  since  sets  of  dimension  zero 
are  rot  permitted,  the  head  is  recogrtz^d  by  the  non-zero  value  of  this 
byte  The  ring  elements  define  rectangular  blocks,  with  the  ring  itself 
representing  the  union  of  :he  rectangles  represented.  There  are  n+2 
half-words  in  each  rectangle  block  the  f’rst  two  being  the  ring  element 
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Table  4.2  Type  Structure  Routi7-,<  s  for  Compiler 


Figure  4.  3  The  set  structur-^  for  n-lupies. 


r  n  i  POINTER  I 
0  78  31 


Figure  4.  4  The  n-tuple  set  head  word. 
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word).  Each  half-word  after  the  first  word  consists  of  two  bytes,  the 
first  of  which  supplies  the  lower  bound  on  the  rectangle  in  one  of  the  di- 
mensiciis,  and  the  second  of  which  supplies  the  upper  bound.  Set  oper¬ 
ators  for  intersection,  union,  and  difference  can  be  readily  designed 
for  this  structure. 

Increments  (the  g-iines  in  events)  will  be  represented  by  the 
same  type  of  head  element  as  for  sets,  showing  the  dimension  and  a 
pointer.  In  this  case,  however,  the  pointer  will  be  to  a  block  of  n 
bytes,  each  giving  the  increment  in  its  correspondir»g  dimension. 

All  blocks  in  these  two  structures  should  be  rounded  to  full  words 
in  their  allocation. 

A  set  of  operators  for  these  two  structures  will  be  required  for 
such  thin  gs  as  intersection,  union,  sum,  etc.  One  special  operator, 
called  back- projection,  is  needed  by  the  compiler.  It  has  two  oper¬ 
ands,  a  set  ard  an  "increme,  t.  "  a.rd  forms  the  structure  for  a  set 
whose  members  are  the  members  of  the  set  operand  m^r-irs  the  value 
of  the  "increment"  operand.  This  set  will  be  represented  by  the  sym- 
bol  v,{S  .  g).  where  S  and  g  are  set  and  increment  oDe’'ands.  respectively. 


5.  THE  COMPILATION  PROCESS 


The  compiler  is  a  program,  called  by  a  single  command  from 
SELMA,  which  trarsforms  the  retwork  structure  into  a  new  structure 
suitable  for  numerical  analysis.  The  new  structure  represents  the 
matrix  U  of  transition  irtens’ties  'see  eq.  2. 1)  in  a  compact  form,  so 
tlut  the  iteration  may  proceed  rapidly  without  extensive  paging  (which 
would  cause  response  to  ihe  "solve”  command  to  be  uncomfortably 
slow'. 

The  compiler  operates  destructively  upon  the  network  structure, 
so  that  a  copy  of  the  network  structure  must  be  saved  prior  to  the 
compilation.  This  furctio'-.  which  calls  upon  the  documentation  phase 
progr«tms,  is  carried  out  by  the  QAS  supervisor  when  it  detecU  a 
’’compile  phase”  command.  Tr  orde”  to  distinguish  between  the  net¬ 
work  structure  and  t^e  worksoace  of  th?  compiler  (which  is  identical 
to  the  retwork  structure  when  rompUitior  begins),  we  will  refer  to 
this  workspace  as  the  working  retwork  area. 

The  compiler  makes  use  of  two  areas,  the  woxking  network 
area  ard  the  trarjit’or.-Ubl^area.  This  latter  cirea  is  the  place  where 
the  new  structure  is  placed  after  the  compilation  is  finished.  It  is, 
obviously,  one  of  the  areas  used  by  the  analysis  procedure  during  the 
analysis  phase  This  a^-ea  is  required  only  at  the  last  stage  of  compila- 
ticn,  ard  is  created  during  the  compile  phase.  •,,If  a  translation  table 
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area  already  exists  at  the  time  a  compile  command  is  given,  its  contents 
will  be  replaced  by  the  new  transition  table.  Hence,  a  documentation 
phase  ’’save  transition  table  area”  command  must  be  given  by  SELMA 
before  the  compile  if  the  old  information  is  to  be  saved. ) 

The  procedure  for  compiling  the  network  involves  a  successive 
absorption  of  the  connections.  Each  connection  is  considered  in  turn, 
and  it  and  the  elements  it  connects  are  replaced  by  a  single  equivalent 
element.  This  process  continues  until  no  connections  remain.  To  illu¬ 
strate,  observe  the  network  of  Figure  5.  la.  The  connections  of  this 
network  will  be  absorbed  in  the  (arbitrary)  order  shown  by  the  ^jncir- 
cled  numbers.  The  compilation  proceeds  as  schematically  indicated  in 
Figure  5,  lb,  c,  d,  e,  fina’ly  resulting  In  a  single  element  having  no  ports, 
which  is  equivalent  to  the  original  network.  This  procedure  provides  a 
systematic  means  for  tracing  the  influence  of  each  au^oevent  upon  the 
entire  network's  behavior.  The  resulting  element  consists  of  a  long 
list  of  autoevents  describing  the  probability  intensity  of  every  c.hange 
in  state  which  is  possible.  The  trinslation  of  this  irfcrmation  Into  the 
form  of  the  transition  table  is  then  a  relatively  sample  operation. 

Figure  5. 2  shows  the  flow  diagram  for  the  compiler.  After  it 
is  invoked  by  a  "compile”  command,  the  compiler  first  determines 
that  no  ports  remain  in  the  network  which  are  not  connected  to  mother 
port,  by  testing  the  ring  whose  head  was  labeled  R  in  the  symbol  table 
of  the  (now)  working  network  structure.  If  this  test  is  successful,  the 
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Figure  5.  2 

The  compiler  flow  diagram. 
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compiler  proceeds  through  a  loop,  successively  absorbing  conrectiors 
until  none  remain. 

The  absorption  of  a  correct'o''  has  been  divided  into  three  dis¬ 
tinct  steps.  In  the  first  step,  the  elements  joined  by  the  connection  to 
be  absorbed  (we  will  call  these  elements  the  associates  of  the  connec¬ 
tion)  are  collected  into  a  single  equivalent  element  having  ports  corres- 
pord'ng  to  each  of  those  of  the  associates.  This  collected  element  then 
replaces  the  associates.  In  the  second  step  the  connection,  which  now 
joirs  c  iy  ports  of  a  single  element,  is  absorbed  to  create  still  another 
elemert.  This  time  the  element  is  equivalent  to  the  collected  element 
with  its  connect ior  ,  and  the  final  replacement  takes  place.  (The  two 
step  operation  eliminates  the  need  to  treat  connections  between  ports  of 
a  s’rgle  element  differently  from  those  between  ports  of  different  ele¬ 
ments.  i  In  the  third  step,  the  sUte  set  is  reduced  to  eliminate  cer^ab 
common  tyijes  of  transient  states 

The  operation  of  the  compMec  ard  its  major  subprograms  is  sum¬ 
marized  in  Table  5.1.  The  rema’^  der  of  this  section  will  d’seuss  the 
subprograms  in  more  detail.  The  subprogram  which  creates  the  ma¬ 
trix  and  state  mapping  structures  will  be  discussed  in  Chapter  6. 

5.  i  C  ollec^tior  of ^Assoejates 

The  operation  called  c  blertior  ol  the  associates  of  a  connect’or 
c  has  a  simple  algebraic  interpretation.  What  it  entails  is  the  crea'ion 
of  a  new  elemert  e’  whose  port  se*  is  the  union  of  the  port  sets  of  .he 
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Table  1  Summ:-.ry  of  Operat’.or.  of  Compiler  aud  Its  Major  Subprograms 


Program  Name  Input  Argumei  t.s  Operations  Performed  Output  Parameters 
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associates,  whose  state  set  is  the  Cartesian  product  of  the  state  sets  ^ 

of  the  associates,  and  whose  events  are  adjusted  to  correspond  to  the 
new  state  set. 

Let  the  set  of  elements  to  be  "collected"  be  E*,  where 

E*  =  {e.:  i€  !♦},  (5.1) 

and  let  the  collected  element  (the  result  of  this  operation)  be  e* 

e*  ^  <r*,p*,P»>  (5.2) 

where 

T*  =  <S*,  H»,Z*  >,  (5.3) 

the  parameter  set  p*  is  null,  and  the  new  port  set  of  e*  is 

i 

P*  =  u  p.  (5.4)  i, 

i  c  I*  ^  i 

I 

Let  I*  =  {ij,  ij, . . . ,  ijj*},  and  select  an  arbitrary  ordering  | 

t 

<  ij)  Xp, . . . ,  i^^  >  on  the  members  of  I*.  The  new  state  set  can  I 

then  be  written  ' 

^  ~  ^  ^  S.  ,  (5.  5) 

where  the  product  x  Ag  of  two  sets  A^,  Ag  of  n-tuples  is  defined 
as  a  set  of  n-tuples  having  a  dimension  which  is  the  sum  of  the  dimen¬ 
sions  of  Aj  and  Ag,  and  for  which  the  members  of  A^  x  Ag  consist  of 
all  concatenations  of  members  of  A^  with  members  of  Ag.  For  exam¬ 
ple, 
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{<1,1>,<2,1>)  X  K3,2>,<3,3>}  =  {<1,1,3, 2>,  <1,1,3, 3>, 

<2,1,3,2>,<2,1,3,3>}  (5.6) 

(Notice  that  this  operator  is  not  precireiy  the  Cartesian  product  operator, 
for  which  the  example  would  give  the  set  {«!,  1>,  <3, 2»,  <<..1, 1>, 

<3,  3»,  «2, 1>  <3, 2»,  «2, 1>  <3,  3»} . ) 

The  autoevent  set  a*  is  the  union  of  autoevents  of  the  elements, 
suitably  modified  to  account  for  the  changed  state  space.  In  particular, 
let  I,  <b.,  g,,  be  an  autoevent  of  an  element  e.  in  £♦  (i.  e.  f  i  L.  , 

k  X.  f.  k  11 

a  a 

i^  •;  I'*').  Then  will  be  represented  by  the  new  autoevent 

1/  =  >,  A?) 

where 


b|*  -  Si  X  ...  X  Si  X  b^  X  Sj  X  ...  X  S,  (5. 8) 

^  h  0-1  *  V+1  n’ 


and 


g^  *  '  <^0, . . . ,  0,  g^ ,  0, . . . ,  0  ^ 

with  g^  concatenated  w  ith  zero-valued  n-tuples  in  a  position  corres¬ 
ponding  to  in  S*.  The  set  S  is 

s»  =  s  Lj,  i£  !•}.  ;5.io; 

To  simplify  notation,  eq.  (5. 8)  will  be  rewritten 

b  *  .  b  (5.11) 

a 

definite  "JCj  "  as  an  operator  meaning  "restrict  the  projection  of  in 
*o 

the  components  corresponding  to  S.  to  the  set  b.".  In  a  similar  veir. 

0  ^ 
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eq.  5. 9  will  be  rewritten 


a 

(5.12) 

defining  "G  "  as  an  operator  meaning  "  extend  g.  into  the  set  S*,  v;ith 

projection  g.  in  the  components  corresponding  to  S.  , 

/y 

and  projection 

U 

zero  in  all  other  components. " 

For  each  port  p^  €  P.,  i  €  !♦.  the  exoevent  function  Z  ♦  is  given, 

similarly  by 

Z*(pj)  =  m  €  M.(pp,  X  e  I*} 

(5.13) 

where 

’  m  m  m  m 

(5.14) 

(5.15) 

(5.16) 

and  where  M.(pp  is  the  index  set,  defined  in  section  4. 3,  of  the  original 
elements.  Each  Z  j*(Pj)  hs  readily  shown  to  be  an  exoevent  set,  satis¬ 
fying  the  constraint  on  given  in  eq.  4. 8,  provided  only  that  ZL(pj) 
was. 

The  collection  routine,  called  by  the  compiler,  makes  the 
changes  in  the  working  network  structure  corresponding  to  the  above 
equations.  Its  input  is  a  pointer  to  a  connection,  and  its  output  is  a 
pointer  to  an  element  e*  which  is  the  "collection''  of  the  connection's 
associates.  The  original  elements  (the  associates)  are  destroyed. 
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This  program,  \vhich  is  flow-charted  in  Fxgjre  5. 3.  has  three  mair 
parts.  The  first  finds  the  associates,  removes  them  from  the  E  rirg. 
and  jo  ns  them  to  a  newly  formed  rirg  E*.  The  second  forms  the  sta*e 
set  S'  of  the  collection  E*  using  fhe  r^tuval  order  of  the  E=^-rirg  as  ♦he 
reverse  order  (right  to  left)  of  tie  prc  jucts  of  the  S.,  i  £  !'».  '.The  or¬ 
der  is  reversed  for  consistency  with  later  operations,  since  the 
set  and  event  sets  and  Z*  will  be  most  easily  Urmf^  by  repeat¬ 
edly  joirlng  successive  components  jmmediately  behind  the  head. )  Tr 
the  course  of  this  activity  the  types  will  be  evaluated  calls  to  the 
’’type  finder"  ( cf.  section  4.  5),  The  third  part  successively  extends 
the  event  sets  and  then  accumulates  both  the  port  sets  and  the  evert 
sets  of  each  element,  destroyir-g  the  elements  as  execution  proceeds 
until  E^  is  empty. 

Figure  5,  4  shows  part  of  the  structure  upon  entry  to  the  coller^o^. 
The  encircled  symbol  represents  a  "bug"  or  pointer  retiired  V  rol- 
lector  program.  Figure  5. 5  shows  that  part  of  the  structure  u.fter  the 
first  two  parts  have  beer  executed  and  the  program  haS  pr^'gressed  to 
the  point  marked  /3  in  the  flow  diagram  ♦Figure  5. 3).  Figure  5.  6  shows 
the  same  part  o!  the  structure  upon  return  from  the  collector.  The 
structure  of  the  state  sets  has  been  omitted  from  the  diagrams  and 
will  be  treated  separately  later,  in  section  5.  6. 

The  first  two  parts  of  this  program  need  no  further  discussicr. 
Pa-'t  three  (point  /3  in  Figure  5.3)  begins  by  taking  the  first  element 


PART  1 


ENTRY 


FIND  ASSOCIATES 


ERROR 


PART  2 


EVALUATE  TYPES  FOR 
ALL  e€E* 

FORM  S^. 


PICK  FIRST  e'€E 
LET  E»-^E»-{e'> 
LET  S'-^S* 

LET  E-^E  U{«'} 
LETp^^  BE  EMPTY 


EXTEND  EVENTS 


PART  3 


EMPTY 


DESTROY  E* 


PICK  First  ecE 


RETURN 


EXTEND  EVENTS' 
i  OFe  > 


LET  P'-*-P'UP 
LET  H*-^a*UH 
LET2'-^2'UE 
DESTROY  e,T 


Figiire  5. 3  Flow  diagram 

for  "called  asso¬ 
ciates"  routine 


ucture  at  li  in  Figure  5 


58 


block  of  E*  and  letting  it  be  the  element  block  which  will  become  the 
collected  element  e*.  It  is  removed  from  E*  and  put  back  in  E.  The 
8=^  which  was  calculated  in  part  two  replaces  the  state  set  of  this  ele¬ 
ment.  The  parameter  list  of  this  element  is  repla.ced  by  an  empty 
list  (the  type  must  henceforth  be  literal).  From  here  on,  following 
programming  convention,  this  element  will  be  referred  to  as  e',  until 
the  end  of  the  operation,  when  it  becomes  the  e*  of  equation  (5. 2).  The 
b^  and  g^  of  the  exoevents  and  autoevents  of  e*  are  extended  to  the  new 

state  space,  and  replaced.  That  is,  b.  is  replaced  by  .  b.  and  g. 

*  a  ^  ^ 

by  8*0.  g  for  each  ^  .  c  H*  and  each  ^  c  Z*.  This  latter  operation 
o  ^  f 

is  done  by  a  separate  subroutine. 

At  this  point  ("a"  in  Figure  5. 3)  each  remaining  element  in  E*  is 
taken,  in  turn,  its  events  similarly  extended  to  8*,  and  its  ports  and 
events  added  appropriately  to  the  ports  and  events  of  e’.  If  P,  H,  Z 
are  the  port  set,  autoevent  set  and  exoevent  function  of  the  current 
element,  then,  since  they  are  each  represented  rings,  they  are 
added  to  P',  H',  Z'  by  insertion 


P’ 

<— 

P’  li  P 

(5.18) 

H  ’ 

<- 

2’  u  H 

(5.19) 

Z' 

<— 

Z'  u  Z 

(5.20) 

where  the  arrow  represents  "replaced  by'*  and  the  union  symbol  here 
means  "insert  the  ring  members  of  the  second  operand  between  the 
head  and  ring  members  of  the  first  operand. "  The  head  of  the  second 
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operand  is  destroyed  in  the  operation,  which  was  called  "ring  concate¬ 
nation"  in  Table  3.  3. 

The  major  subprograms  of  the  collect  routine  are  summarized 
in  Table  5. 2. 

5. 2  The  Meaning  of  Connections 

The  exocvents  at  the  ports  of  an  element  represent  the  conditi'  ni 
under  which  tasks  can  be  removed  from  (or  inserted  into)  the  <;lement 
via  the  port,  and  the  consequence  to  the  element's  state  of  that  removal 
(or  insertion).  The  function  of  a  connection  is  to  indicate  a  relationship 
between  the  exoevents  of  a  set  of  ports.  It  indicates  the  conditions  under 
which  a.  task  Ls  to  be  transferred  between  a  pair  of  ports.  For  example, 
when  a  pair  of  ports  is  joined  by  a  simple  connection,  a  transfer  of  a 
task  occurs  immediately  whenever  the  "outpul"  port  has  a  task  which  can 
be  removed  and  the  "input"  port  simultaneously  can  accept  a  task.  The 
exoevents  at  the  two  ports  are  thought  to  "occur"  s;x)ntaneously  wherever 
these  conditions  are  met.  For  each  type  of  connect'on  the  rules  deter¬ 
mining  the  occurrence  of  the  exoevents  will  differ,  and  the  meaning  of 
connections  is  conveyed  by  these  rules. 

Since  we  will  be  operating  onJly  on  collected  elements,  the  only 
time  we  are  concerned  with  the  meaning  of  a  connection  will  be  when 
it  is  a  connection  between  ports  belonging  to  a  single  element.  The  def¬ 
initions  that  follow  will  assume  that  condition. 


' 


i  Mi 


Table  5.  2  Summary  •'-'Associate  Collection  Routine 
(continued) 
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The  information  describing  the  effect  of  such  a  connection  upon  its 
associated  element  will  be  described  by  a  set  of  objects  called  spontane¬ 
ous  events.  These  events  describe  the  conditions  for,  and  immediate 
consequences  of,  transfer  of  tasks  through  the  connection.  These 
events  can  be  thought  to  occur  "spontaneously”  whenever  the  state  of 
the  element  becomes  one  of  a  set  of  forbidden  states.  For  example,  in 
the  simple  connection,  a  forbidden  state  would  be  one  in  which  one  port 
was  "offering"  a  task  and  the  other  could  "accept"  one.  The  consequence 
would  immediately  be  a  further  change  in  state. 

Let  the  spontaneous  event  set  of  a  connection  Cj^  joining  ports 
of  the  element  e^  be 

=  {9^:htH^).  (5.21) 

where 

b.  is  a  subset  of  S. 
h  1 

gj^  is  a  constant  n-tuple  such  that 
for  each  x  €  x  gjj  is  i  i  S. 

7^  is  a  probability 
is  an  index  set. 

The  set  bj^  is  the  forbidden  state  set,  gj^  is  the  increment  of  the 
event,  and  is  the  probability  given  that  the  state  is  in  that  it  will 
jump  by  an  amount  gj^.  As  fur  exoevent  sets,  the  probabilities  in  a 
spontaneous  event  set  must  satisfy  the  restriction  that,  for  each  x  c  S. 


/ 
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2  TT.  -=Oor.l.  (5.23) 

h:  " 

Xfbh 


The  spontaneous  events  are  formed  from  the  information  con¬ 
veyed  in  the  exoevent  sets  of  the  connected  ports  in  a  manner  which 
is  distinct  to  the  '’type"  of  the  connection.  Let  a  connection  Cj^  be  of 
simple  type.  Let  the  ports  joined  by  Cj^  be  -  {p.  ,  p.  }  and  let 
those  both  be  ports  of  the  element  e.,  so  that  P.  C  P..  ilecall  that 
the  exoevent  sets  of  p.  and  p.  are  given  by 

h  b 


-  “l  ^ 

(5. 24) 

“  “a  e  Mj(Pj^)} 

{5.  25) 

m* 

1 

“l  “l 

(5.  26) 

( 

“2 

^  b  ,  g  ,  >  . 

“2  “2  “2 

f.5. 27} 

The  spontaneous  event  set  6^  for  this  simple  connection  is  given  by 

ejj  {9^:  h  ^  Mj(Pj  )  ;<  Mj(Pj  )}  (5. 28) 

X  2 

where 


This  indicates  that  spontaneous  events  occur  whenever  the  state  is  in 
both  endocondition  sets,  and  that  the  increment  is  the  sum  of  that  due 
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to  removing  a  task  from  one  port  and  that  due  to  inserting  a  task  into 
the  other  port.  The  set  of  states  for  which  no  spontaneous  event  occurs 


b,  =  S.  -  u  (b  n  b  ). 
m^,m2  1  2 


(5.  30) 


For  convenience  in  the  programming,  we  will  find  an  augmented 
spontaneous  event  set  to  be  useful.  It  is  defined,  for  the  simple  con¬ 


nection,  as 


®k  =  ®k  ^  0,  1  >}, 


(5.  31) 


where  the  additional  event  represents  no  change  for  states  which  are  not 
forbidden.  In  this  case,  the  probabilities  add  to  1  for  every  member  of 
Sj.  It  is  readily  shown  that  eq.  5. 23  is  satisfied  for  the  spontaneous 
event  set  ej^^. 

The  other  connection  types  will  be  discussed  below.  In  what  fol¬ 
lows  the  element  whose  ports  are  being  connected  will  be  called  e.,  and 
the  exoevent  set  of  any  one  of  its  ports  p.  will  be  assumed  to  be  written 


(S.  32) 


a  a  a  a 


(5. 33) 


The  probability  summation  property  will  hold  for  all  spontaneous  event 
sets  defined. 
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Considering  a  connection  Cj^  of  overflow  type,  let  tiie  output  port 
joined  by  c.  be  p.  ,  the  preferred  input  port  be  p.  ,  and  the  overflow 
input  be  p.  ,  so  that  P.  •  {p.  ,p.  ,p.  }  C  P,.  The  spontaneous  events 

Jq  -^2  '  ■ 

consist  of  two  classes*  those  which  represent  a  preferred  task  flow, 
and  those  which  represent  an  overflow.  The  first  class  occurs  when¬ 
ever  the  state  is  in  both  b_  and  b_  ,  for  each  and  m,  in  M(p.  ) 

“O  “l  0  1  Jq 

and  M{p.  )  respectively,  since  that  is  the  only  condition  where  the  pre- 

ferred  flow  is  possible.  The  second  class  occurs  whenever  the  state 

is  in  both  b_  and  b_  ,  but  not  in  anv  of  the  b  ,  m,  s  M(p.  )  (i.  e., 

preferred  flow  is  not  possible),  for  each  and  m«  in  M(p.  )  and  M(p.  ), 

w  ^  Jq  Jj' 

Algebraically, 

^  M{Pj  )  X  M(p.  )}  u  h  t:  M(p.  )  X  M(p.  )}  (5.  34) 


“0  “0 


(5.  35' 


t>  _  ^  0  n  b  -  (  0  b 

“0  “2  m,  -M(p.  '“l 

1 


“1  “0  “2  “2 


(5. 36) 


for  M(p.  ),  m.  t  M(p.  ),  m,  ;  MCp.  ). 

Jq  ^  ^2 


The  ro-event  state  set 


b.  S.  -  u  b. , 
^  h  :  H  ” 


(5. 37) 


where  H  (M(p.  )  x  M(p.  ))  u  (M(p.  ;  Mlp.  )),  and  the  augmented 

Jo  n  ’0  ’2 

spontaneous  event  set 


©t  -  ®ic  ^  {<b,,0,l  >}. 


{5. 38) 


t 
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The  merge  type  connection,  joining  an  input  port  p^  of  e^,  a  (pre¬ 
ferred)  output  port  pj ,  and  (secondary)  output  port  Pg  (oae  must,  in  de¬ 
fining  any  merging-type  elements  resolve  the  competition  among  inputs, 
and  that  is  done  here  on  a  purely  priority  basis)  has  precisely  the 
same  spontaneous  events  as  the  overflow -type  element.  Thus  they  are 
indistinguishable,  and  will  henceforth  be  given  the  same  type-name,  the 
"overflow -mergft  " 

The  random  branch  is  defined,  for  the  purposes  of  this  report, 
as  a  branch  for  which  flow  can  occur  only  if  a  task  can  be  accepted  at 
any  connected  input  port  of  the  element,  and  if  a  task  is  also  available 
from  the  connected  output  port.  Let  the  output  port  of  the  element  be 
Pq  and  the  input  ports  be  p^, . . .,  p^,  where  N  is  the  number  of  branches 
in  the  connection.  When  this  condition  is  met  a  task  is  transferred 
from  Pq  to  one  of  the  ports,  p.  say,  with  probability  y  ,  o=l,  2, , . . ,  N, 
where  (N,'/^, . . .  ,7^^)  is  the  parameter  list  of  the  connection.  Alge¬ 
braically, 


®k  ®  ^^^0^  X  ...  X  M(Jj^)  X  {1, . . . ,  n}} 


“0 


‘N 


(5.  39) 

jTT  •••7T 

(5,  40) 

for  all  m  €  M(p.  o  =  1, . . . ,  N.  As  usual 


(5.41) 


and 
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\  -  QfcU  {<bj^,0,l  >}.  (5.42) 

Provided  that  the  property 

f  -  1  (5.«) 


holds,  the  probability  summation  property  (eq.  5. 23)  will  also  hold  for 


5. 3  The  Absort^ion  of  Comiections 

The  operation  of  connection  absorption  removes  the  ports  con¬ 
nected  from  the  associated  element,  prunes  the  exoevent  sets  at  the 
removed  ports  from  the  exoevent  function,  forms  the  spontaneous 
event  set,  and  then  consolidates  the  autoevents,  remaining  exoevents, 
and  spontaneous  events,  to  produce  a  new  element.  The  operation  of 
this  subroutine  of  the  compiler  is  illustrated  schematically  in  Fi^re 
5. 7,  and  its  major  parts  are  summarized  in  Table  5. 3. 

The  first  operation  in  the  program  determines  the  element  joined, 
and  alters  the  structure  in  preparation  for  the  formation  of  the  spon¬ 
taneous  events.  The  ports  joined  by  the  connection  are  removed,  and 
appropriate  alteration  of  the  exoevent  function  is  carried  out,  with  the 
exoevent  sets  of  the  removed  ports  shifted  to  the  connection  structure. 
This  operation  is  clearly  illustrated  in  Figures  5. 8  and  5. 9,  which  show 
fragments  of  the  network  structure  before  and  after  the  first  block  of 
Figure  5. 7  is  executed  In  the  course  of  finding  the  element  associated 


oo 


Table  5.  3  Summary —Absorb  Connection  Routine 
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Figure  5.  8 

Fragment  of  working  network  structure  at  beginning 
of  connection  absorption. 
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with  the  connection,  this  block  will  also  dlsco?ei  if,  by  failure  of  the 
element  collection  program,  the  number  of  elements  connected  is  not 
unity.  This  is  a  fatal  error. 

At  the  conclusion  of  these  operations,  an  appropriate  Spontaneous 
Event  Routine  is  called,  depending  on  the  connection  type.  Since  the 
spontaneous  events  will  not  be  altered  during  the  period  of  their  exis¬ 
tence,  and  since  the  operation  of  "getting”  and  "freeing"  small  blocks  is 
relatively  time-consuming,  the  Spontaneous  Event  Routines  should  create 
the  entire  6 -structure  so  that  it  is  physically  contiguous.  If  this  is 
done,  only  one  fairly  large  block  need  be  gotten  from  the  free  core 
manager,  and  it  will  be  kept  until  the  connection  absorption  is  complete, 
at  which  time  the  entire  structure  can  be  destroyed  by  a  single  "free" 
command  to  the  manager  (see  "Destroy  ia  Figure  5. 7). 

The  spontaneous  event  routine  (a  prototype  of  whfeh  is  shown  in 
Figure  5. 10)  constructs  a  list  structure  which  correspomu.  to  the  set 
described  in  the  previous  section  in  its  requested  core  space.  The 
events  will  be  generated  one  by  one  and  added  to  a  ring,  while  the 


state  set  and  no-event  set  bj^  will  be  pruned  with  the  formation  of  each 


event.  Thus,  at  the  conclusion  of  the  formation  of  the  last  spontaneous 


event,  the  state  set  has  had  all  forbidden  states  removed  from  it,  and 


the  augmented  spontaneous  event  set  6^^  has  been  formed.  Whenever  a 
forbidden  state  set  b|^  is  found  to  be  empty,  the  corresponding  spontaneous 
event  will  be  simply  omitted.  In  this  routine,  also,  the  parameters 


Figur  ,*  5. 10 

Prototype  fiow  diagram  for  spontaneous  event  routines. 
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are  checked,  if  this  is  a  random  branch,  to  see  that  the  probabilities  (y^ 
in  section  5.  2)  sum  to  unity.  The  exoevent  sets  for  the  connected  ports 
are  being  destroyed  after  the  spontaneous  events  have  been  formed  ,  and 
the  structure  finally  appears  as  illustrated  in  Figure  5. 11. 

5. 4  Consolidation  of  Events 

The  operation  of  event  consolidation  replaces  each  event  (either  exo¬ 
event  or  autoevent)  in  the  element  e.  associated  with  the  connection  by  a 
new  event,  or  events,  which  includes  the  effect  of  the  spontaneous  events. 
The  action  to  be  accounted  for  in  this  operation  is  that  produced  when  an 
event  of  the  element  causes  a  spontaneous  event  to  occur. 

Such  is  the  case,  for  example,  when  a  task  completion  in  a  server 
causes  a  transfer  of  the  completed  task  through  the  connection  joining  its 
output  port.  The  effect  of  the  autoevent  is  to  change  the  state.  If  the 
state  is  thereby  changed  to  a  forbidden  state,  the  spontaneous  event 
(task  transfer)  occurs,  changing  the  state  again.  In  this  report,  we 
assume  that  only  one  task  may  be  required  to  pass  thrcjgb  a  connection 
upon  occurrence  of  any  event,  so  that  the  second  state  charge  cannot 
itself  result  in  another  spontaneous  event  at  the  same  corjiection. 

A  similar  argument  applies  to  exoevents,  as,  for  example,  the 
consequences  of  a  task  being  entered  into  the  input  of  a  queue.  This 
will  result  in  an  increase  in  state  which  may  cause  a  trar^sfer  of  a  task 
through  the  connection  joinmg  its  output.  That  is,  the  state  is  chained 


77 


b  ,  g 
m. 

J  j 


>  > 
m. 

J 


(5.  51) 


becomes 


where 


«*m.h  “  <  *’m.  '■  \  <V’  8m.  8h>  %.  ’'h  >  '  '5' 

1  1  j  j 


(The  probability  property,  eq.  4. 18,  for  exo-event  sets  is  readily 
shown  to  be  preserved. ) 


be  empty.  In  such  a  case,  the  event  conveys  no  us^ul  information  and 
may  be  omitted  from  the  corresponding  set  of  events. 

The  program  which  will  accomplish  this  transformation  is  the 
event  consolidator,  described  by  the  flow  diagram  of  Figure  5. 12. 

This  program  calls  on  the  "event  set  consolidator"  which  is  described 
the  diagram  of  Figure  5. 13.  The  symbol  $  is  used  in  this  diagram 
to  represent  an  arbitrary  event  set 

$  d  €  D}  (5.  53) 

♦d  ’  ^'’d’8d’''d!  • 


This  program  proceeds  around  the  ring,  inserting  the  new  events 
behind  it  as  it  proceeds.  Each  time  a  has  been  treated  completely, 
it  is  destroyed  and  the  program  proceeds  to  the  next  until  it  comes 


tYES 

RETURN^ 

Figure  5. 12 

Flow  diagram:  Event  consolidator. 
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back  to  the  head  of  the  ring.  If  either  set  is  empty,  the  result  is  an 
empty  $♦. 

5.  5  Trimmit^  of  the  State  Set 

The  consolidation  procedure  described  above  can,  under  certain 
circumstances,  produce  surprisingly  large  state  sets.  This  will  be  par¬ 
ticularly  troublesome  when  a  random  branch  was  involved  in  the  compila¬ 
tion  at  some  stage.  This  is  an  manifestation  of  a  general  problem  facing 
this  translation  which  cannot  be  economically  solved:  the  problem  of 
transient  states  and  reducible  Markov  chain  models. 

The  points  in  the  rectangular  "state  set"  Sj^  of  the  (finite)  Markov 
chain  of  the  entire  network  fall  basically  into  three  categories 

1)  Forbidden  states 

2)  Recurrent  states 

3)  Transient  states. 

The  forbidden  states  are  points  of  which  result  immediately  in  spon¬ 
taneous  events,  which,  in  turn,  take  the  process  to  a  point  which  is  not  a 
forbidden  state.  The  forbidden  states  are  automatically  pruned  from  the 
state  set  upon  veach  consolidation,  and  hence  are  no  problem. 

The  present  difficulty  revolves  around  the  presence  ol  transient 
states  in  the  modeL  Transient  states  are  legitimate  states  of  the  Maikov 
model.  However,  since  they  have  an  equilibrium  probability  of  zero, 
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and  since  information  bearing  upon  them  ordinarily  has  no  bearing  on  the 
calculation  of  equilibrium  probabilities  of  other  states,  they  are  "excess 
baggage"  and  will  waste  valuable  storage  in  the  various  data  structures 
of  QAS. 

The  ideal  procedure  would  be  to  identify  all  transient  states  of  con> 
solidated  elements  as  they  are  created,  and  to  remove  them  from  the 
state  set.  Unfortunately,  there  appears  to  be  no  systematic  procedure 
for  identifying  all  of  the  transient  states  without  performing  a  calculation 
very  similar  to  evaluation  of  all  equilibrium  probabilities,  which  defeats 
our  purpose.  (We  want  the  state  set  reduced  before  we  calculate  equili¬ 
brium  probabilities. ) 

The  example  of  the  random  branch  is  a  good  one.  Let  a  diagram 
contain  a  fragment  like  that  of  Figure  5. 14.  The  states  for  which  more 


Figure  5. 14 

A  network  fragment  containing  a  random  branch. 
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than  one  server  is  busy  cannot  be  entered  from  any  state  for  which  less 
than  two  servers  arc  busy,  since  the  random  branch  is  blocked  when  one 
server  is  busy,  and  no  "tasks"  will  be  transferred.  However,  if  the  sys¬ 
tem  starts  with  more  than  one  server  busy  it  will,  upon  exit  of  tasks 
from  service,  enter  states  with  less  than  two  servers  busy.  Thus,  there 
is  a  hiirly  large  number  of  transient  states  (37  recurrent  states,  120  tran¬ 
sient  states,  4  forbidden  states).  As  more  elements  are  consolidated  with 
this  fragment,  the  ratio  of  transient  to  recurrent  states  will  remain  much 
the  same,  at  best  decreasing  to  2:1.  The  states  of  the  example  were  easy 
to  identify  through  physical  insight,  a  weapon  not  available  to  the  compiler. 

A  general  solution  to  the  problem  of  eliminating  a  transient  states 
is  not  available.  A  more  specific  "fix"  must  be  one  which  is  guaranteed 
valid  for  every  network  and  every  definable  element  and  connection.  In 
other  words,  it  should  not  unnecessarily  restrict  the  already  restricted 
class  of  models  which  can  be  treated  by  QAS.  Especially,  it  should  not 
result  in  a  restriction  on  inputs  to  the  compiler  which  cannot  be  enforced 
syntactically  before  compilation  is  attempted. 

Fortunately,  there  is  a  class  of  transient  states  which  is  relatively 
easily  recognized  and  which  includes  those  of  the  example.  The  recom¬ 
mended  solution  will  eliminate  this  particular  kind  of  transient  state, 
and  result  in  perfect  pruning  of  states  for  the  random  branch.  Figure 
5. 1 5  illustrates  various  cases  of  communication  among  transient  (shown 
by  circles)  and  recurrent  states  (shown  by  x  marks)  for  a  total  of  three 


. . 
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states. 


.  ■  ji)  nci-  on  ~  c.r-  nt  s 

S'  .<wn  ...  ••  n-:  o.  ) 


The  procedure  recommended  is  that,  subsequent  to  the  absorp¬ 
tion  of  each  connection,  the  states  which  are  not  "target  states"  of  any 
event  be  successively  trimmed  from  the  state  set  until  no  more  such 
states  can  be  found.  Since  case  (a)  of  Figure  5. 15  is  representative 
of  what  happens  in  the  random  branch  example,  the  need  for  success¬ 
ive  trimming  should  be  clear.  The  target  states  of  an  event  6^  (either 
auto-  or  exo-l  are  found  by  the  expression 

n  (b.) 

-«d  ^ 


using  the  "oack  projection"  operator  j;  previously  defined. 
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Figures  5. 16  and  5. 17  illustrate  the  programs  required,  and  show 
how  the  removal  of  previously  located  transient  states  can  be  efficiently 
accomplished  while  new  transient  states  are  being  sought.  The  set 
represents  a  set  of  known  transient  states,  while  the  set  T2  is  a  set  of 
states  formed  by  starting  with  the  state  set  and  repeatedly  deleting  tar¬ 
get  state  of  events  as  they  are  sequenced  through.  Once  Tg  is  found  to 
be  empty,  we  know  that  no  further  transient  states  will  be  found  by  this 
method  However,  the  sequence  must  be  continued  if  is  not  empty, 
so  that  all  previously  located  transient  states  are  removed  from  all 
events.  Figure  5. 17  is  not  particularly  efficient  as  shown,  since  cal¬ 
culations  of  set-differences,  and  evaluations  of  ,7  are  performed  even 
when  some  components  (T^  or  Tg)  are  known  to  be  empty.  Suitable 
switch  setting  or  branching  can  make  this  program  highly  efficient  since 
the  most  probable  situation  is  expected  to  be  the  case  with  empty, 
and  since  Tg  is  expected  to  be  frequently  empty  after  only  a  few  events 
are  examined  Figure  5. 18  shows  one  possible  configuration  which  will 
take  ad\'antage  of  such  information. 

5.  6  The  Result  of  Collection,  Absorption  and  Trimming 

After  repeated  collection  of  associates  and  absorption  of  connec¬ 
tions,  the  working  network  structure  will  eventually  lack  remaining 
connections.  At  that  point,  the  network  structure  will  ordinarily  look 
like  Figure  5, 19 ,  with  a  single  element  having  no  ports  or  exoevents, 


TEST  AND  \ 
TRIM  EVENTS 
\0F  £|(p) 


RETURN 


Figure  5. 16 

State  sei  trimmer  flow  diagram. 
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entry) 


Figure  5. 17 

"Test  and  trim  events”  flow  diagram. 


PICK  FIRST 


Figure  5. 18 


A  more  efficient  ’’test  and  trim” 


program. 
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and  a  large  set  of  autoevents.  The  information  in  this  set  of  autoevents 
describes  the  information  in  the  matrix  of  transition  intensities  for  the 
netivork,  which  must  be  prepared  next. 

It  is,  however,  possible  that  the  network  will  consist  of  more  than 
one  element,  with  all  the  elements  isolated  from  one  another.  There 
can  still  be  no  ports,  since  there  are  no  connections,  and  there  were 
originally  no  unconnected  ports.  This  multi-element  condition  is  the 
suit  of  the  original  network  consisting  of  two  or  more  isolated  parts. 
This  im’^lies  that  the  model  described  was,  in  fact,  several  models, 
each  of  which  has  been  compiled  here.  At  present,  this  situation  can 
be  treated  as  an  error  condition.  Nevertheless,  it  offers  interesting 
potentialities  for  more  flexible  use  of  the  system,  and  future  implemen¬ 
tation  should  probably  allow  for  it.  In  that  case^  procedures  would  be 
needed  for  identifying  which  isolated  network  was  to  be  solved,  and  for 
what  variables.  For  the  time  being,  however,  it  is  easier  to  give  an 
error  return  if  the  result  of  absorbing  all  connections  is  a  network  with 


more  than  one  element. 


6.  THE  STATE  AND  TRANSITION-MATRIX  STRUCTURES 


In  the  numerical  calculation  of  equilibrium  probabilities,  it  is  neces¬ 
sary  to  repeatedly  multiply  the  "matrix”  of  transition  intensities  by  a 
probability  "vector. "  The  autoevent  set  resulting  from  the  repeated 
absorption  essentially  describes  the  matrix  of  transition  intensities:  the 

set  b,  describes  a  set  of  columns  of  the  matrix  which  contain  the  inten¬ 
se 

sity  while  describes  how  far  off-diagonal  the  intensity  is  in 
each  of  these  columns — hence  defining  the  row  in  which  each  can  be 
found  The  state  set  resulting  from  the  repeated  absorption  describes 
the  index  set  of  the  probability  vectors.  (To  each  state  there  must 
correspond  a  unique  probability  in  the  probability  vector). 

Because  of  the  multidimensional  description  of  states,  however, 
the  structure  that  results  from  the  repeated  absorption  is  unsuited  to 
rapid  execution  of  the  matrix-vector  multiplication.  To  use  this  struc¬ 
ture  in  its  original  form  one  must  either  spread  the  probability  vector 
over  a  prohibitively  large  area  of  memory  so  that  a  very  simple  formula 
can  be  used  to  find  a  probability  when  a  state  is  given,  or  one  must 
spend  enormous  amounts  of  time  in  the  process  of  state  mapping  re¬ 
peatedly  upon  each  multiplication.  The  large  number  of  dimensions 
involved  makes  it  likely  that  a  simple  linear  mapping  of  a  large  rec¬ 
tangular  prism  containing  all  the  states  may  actually  use  hundreds  or 


even  thousands  of  times  as  many  locations  as  there  are  states.  Even 
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so,  that  mapping  would  have  to  be  performed  twice  with  each  scalar 
multiplication  in  the  matrix  multiplication,  and  would  be  dominant  in 
determming  the  time  required  for  the  operatior.  Any  more  efficient 
mapping  would  require  much  more  time  and  be  completely  unreason¬ 
able.  Such  time  demands  would  defeat  the  primary  purpose  of  this 
project,  which  is  to  provide  sufficiently  rapid  (and  inexpensive)  res¬ 
ponse  to  enable  the  user  to  experiment  with  models  on  the  computer 
in  a  symbiotic  manner. 

As  a  consequence,  the  only  alternative  is  to  prepare  a  new  struc¬ 
ture  which  is  well  suited  to  the  matrix-vector  multiplication-operation 
operation,  and  permits  rapid  execution  with  minimal  storage  require¬ 
ments. 

The  expense  incurred  is  the  addition  of  another  stage  of  transla¬ 
tion  between  the  consolidated  network  structure  and  this  new  structure. 
This  operation  is  a  part  of  the  compiler  and  was  shown  in  the  flow  dia¬ 
gram  in  Figure  5. 2.  The  objective  is  to  allow  the  matrix- vector  mul¬ 
tiplication  to  proceed  with  such  computational  efficiency  that  the  time 
taken  is  only  ten  to  twenty  percent  greater  thar.  that  taken  by  the  scalar 
multiplications  and  additions  which  are  inherent  in  the  matrix  operation. 
This  section  will  describe  first  the  state  mapping  function  and  the 
structure  which  supports  it.  It  will  then  describe  the  structure  to  be 
used  to  represent  the  matrix.  Section  6. 3  will  finally  describe  the 
procedure  used  to  create  these  structures. 
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6. 1  The  State  Mapping  Function 

To  obtain  an  efficient  use  of  memory  for  the  probability  vectors, 
the  state  set  should  be  mapped  onto  a  consecutive  set  of  integers. 

This  has  the  advantage  of  minimizing  the  number  of  pages  to  be  used 
for  the  probability  vectors,  and  hence  decreasing  the  elapsed  time 
(rather  than  epu  time)  of  solution  under  dynamic  paging  in  the  computer 
The  state  set  has  been  represented  as  a  union  of  rectangles  (in  n  diraen- 
siors).  Let  the  state  set  to  be  mapped  be  S,  and  let 

S  .  Sj  u  u  ...  uSj^  (6.1) 


where  S  is  a  rectangle, 


s  K  . v  }  ,  (v  ....,V  } 

T1  *12  ^21  ^22 


rl  r2 


(6.  2) 


Each  V  is  an  integer,  c^rrespordTg  to  appropriate  values  in  the  n- 
tuple  structure  and  n  is  the  d*mers'or  of  the  set,  for  1,  2, ... ,  Ng. 
The  mapping  will  associate  a  unique  integer  with  every  state  by 


providVg  a  simple  expansion  of  each  rectangle.  Let  6  be  the  car- 

% 


dinality  of  the  union  of  the  rectangles  Sj. .  . .  - 


'  i-l  1  '-1 

^  ^  *-  S;  ^  (S.) 

0  3  1  ' 


(6.  3) 


(6.  4) 


and  let 


0. 


Then  6^+1  represents  the  state  number  corresponding  to  the  lowest 
numbered  state  in  The  numbering  of  the  remaining  states  within 
is  according  to  the  usual  rectangular  mapping  used  in  most  programming 
systems.  Let  6  be  the  product  of  cardinalities  of  the  previous  i-1  di¬ 
mentions  of  S  ,  so  that 


6  =  1 


5  -  5  (v  -V  +1),  1-1,...,  n. 


(6.  5) 

(6.  6) 


“'il 


Then,  accordingly,  a  vector  x  =  <  x^, . . . ,  >  which  is  a  state  in 

will  be  mapped  to  the  integer 


X  r-  6  (x,  -  V  )6  + 


. . .  -j-  (x  -  V  )6 

'O  *  'll  "I  "  *nl  "n 


(6.7) 


In  this  manner,  each  of  the  states  in  S  has  a  distinct  mapping  x 

and  the  states  in  S  map  to  a  set  of  consecutive  integers,  {l, . . . ,  6  q}. 

S 

Complete  information  required  to  us-"’  the  state  mapping,  Eq.  <6.  7), 
is  giver  by  the  tv'o-dimer.sional  array  <  5  ,  :  v-1,, . . ,  N„,  '-0, 1, . . . ,  n  >. 

3 

The  state  mapping  structure  produced  by  the  compiler  is  simply  a  MAD 
(or  FORTRAN) -type  two-dimensional  array  of  halfword  integers,  as¬ 


suming  (reasonably)  that  the  number  of  states  will  always  be  less  than 


65,  536. 


Operators  required  to  construct  this  structure  include  one  (a 
’’state  structure  generator” .  which  prepares  the  structure  when  the 
state  set  in  n-tuple  set  form  is  given,  and  another  (a  ’’state  mapper”) 
which  uses  this  matrix  and  the  n-tuple  form  of  S  to  obtain  the  mapped 
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value  of  a  given  state  according  to  Eq.  6.  7  . 

6.2  The  Transition  Matrix  Description 

The  matrix  structure  must  describe  the  autoevents  H  in  such  a 
form  that  no  state  mapping  need  be  performed  during  matrix-vector 
multiplication.  It  should  permit  rapid  execution  and  efficient  storage 
for  the  common  forms  encountered. 

The  requirements  seem  to  demand  that  the  structure  be  in  some 
type  of  outline  form  whereby  events  which  occur  for  rectangular  sub¬ 
sets  of  the  state  sets  can  be  compactly  represented  and  evaluated 
through  direct  indexing  on  the  data  of  the  structure.  One  such  struc¬ 
ture  was  proposed  in  [3]  for  a  somewhat  more  general  class  of  prob¬ 
lems,  and  that  reference  gives  srme  useful  ideas,  A  more  special¬ 
ized  structure  which  appears  to  be  more  efficient  will  be  incompletely 
presorted  here. 

The  structure  is  illustrated  in  Figure  5. 1.  A  value-lire  or 
in  the  figure)  indicates  the  value  of  co’lecticr  of  idertically  valued 
matrix  elements.  A  blcc^-line  g  'r  the  figure)  indicates  the  first 
row  ’rdex  of  a  subset  ;f.alled  a  "block”;  of  the  collection  of  iden¬ 
tical  matrix  elements,  and  a  number  (gi  which,  when  added  to  a  column 
index,  gives  a  corresponding  row  irdex  fo’'  each  member  of  the  block 
of  matrix  elements.  The  members  of  the  block  are  identified  by  a 
series  of  indexing-lines  t/3, 6  in  *he  figure),  which  indicate  increments 
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Figure  6. 1  Illustrating  a  matrix  structure. 

(6  )  to  be  applied  to  the  column  index  to  get  another  column  index, 
and  counts  {li)  of  the  number  of  times  the  increment  is  to  be  applied. 
Each  line  is  a  word,  the  /^q’s.  6 's.  and  g’s  are  nalfwords.  the  ,;'s  are 
bytes,  and  the  remaining  bytes  are  coded  with  sufficient  information 
so  that  line  types  can  be  recognized. 

This  structure,  when  the  details  are  worked  out.  should  be 
capable  of  describing  any  matrix  at  least  as  efficiently  as  by  a  list 
of  triples  (value,  row  index,  column  index),  and  much  more  efficiently 
when  element  values  are  repeated  diagonally  in  a  regular  pattern,  as 
will  be  the  case  for  QAS  transition  intensity  matrices  using  the  above- 
described  state  mapping  procedure.  Because  tne  information  con¬ 
tained  in  the  indexing-lines  is  exactly  the  information  required  for 
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quickly  loading  ano  testing  "index-registers.  "  rapio  execution  of  tne 
matrix-vector  multiplication  should  also  be  assured. 


To  further  illustrate,  consider  the  matrix  U: 


where  only  two  sets  of  elements  of  the  matrix  are  shown.  These 
elements  would  be  described  in  the  matrix  structure  as  in  Figure  6.  2. 
The  first  block  is  the  set  of  values  ”//’  in  the  fourth,  fifth,  seventn. 
eighth,  tenth,  eleventn,  thirteenth,  and  fourteenth  rows.  Thus  tne 
/^Q  index  is  four.  Each  element  of  this  block  is  in  a  column  whose 
index  is  two  less  than  its  row.  Thus  the  g  for  this  block  is  -2.  The 
value  repeats  twice  with  unit  spacing,  and  that  pair  repeats  four  times 
with  spacing  of  three.  Two  other  blocks  consisting  of  single  elements 
(i.  e.  b-lines  only)  complete  the  treatment  for  /x.  For  there  is  a 
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Figure  6.  2  The  structure  for  the  matrix  in  eq.  (6.  9). 
single  block  starting  in  the  first  row  (/5q=1-  i=5.  etc. ). 

Multiplication  of  this  matrix  by  a  row  vector  proceeds  by  con¬ 
sidering  each  block  of  tlie  matrix  structure  in  succession.  The  value 
pL  is  placed  in  a  register  and  index  registers  are  set  up  for  nested 
loops,  with  depth  of  nesting  equal  to  the  number  of  i-lines  in  the  block. 
The  gives  the  initial  displacement  of  index  of  the  appropriate  ele¬ 
ment  on  the  vector  to  be  multiplied  by  pt.  and  g  indicates  the  increment 
to  be  used  to  locate  the  i.ndex  of  the  element  of  the  result  vector  into 
which  the  product  is  to  be  accumulated. 

To  be  more  concise,  if  y  and  x  are  vectors  ana 
y  =  xU, 

and  if  a  vector  y  is  initially  zero,  then  yii-Q*  g  +  ij  I2  ^  ^ 
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replaced  by 

repeatedly  for  each  -  0, 6  25  j, . . . ,  /3j6  ^ ;  ig-^O,  6  g,  26  g?  • . .  >  ^2^2’ 

. , .  ;  i  0, 6^  26  , . . . , /3  6 where  q  is  the  number  of  i-lines  in  the 

q  q  q  q  Q 

block.  This  replacement  is  repeated  for  each  block  of  the  matrix  struc¬ 
ture.  The  result  is  the  vector  y. 

By  suitable  use  of  registers »  this  operation  should  not  require 
significantly  more  time  for  execution  than  is  taken  by  the  scalar  mul¬ 
tiplications  and  additions  alone. 

6. 3  Creation  of  Transition- Matrix  and  State  Mapping  Structures 
The  proposed  outline  structure  for  matri :  representation  is 
well  suited  to  the  storage  of  transition  intercity  matrices  because  it 
is  efficient  for  both  isolated  matrix  elements  and  for  nested,  diag¬ 
onally  repetitive  elements.  These  two  rases  occur  with  great  fre¬ 
quency  in  transition  intensity  matrices  of  Markovian  queueing  net¬ 
works. 

The  final  subprogram  of  the  compiler  has  the  task  of  providing 
the  state  mapping  structure  and  the  transition  matrix  structure. 

These  are  created  in  two  separate  areas:  the  state  area,  and  the 
matrix  area.  These  areas  are  created  by  suitable  calls  to  the  QAS 


supervisor. 
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The  state  and  matrix  constructor  (a  subprogram  of  the  compiler) 
contains  two  subprograms.  The  first  creates  the  state  area  and  its 
structures.  The  second  converts  an  autoevent  into  a  series  of  lines  in 
the  matrix  structure.  The  state  and  matrix  constructor  calls  the  lat¬ 
ter  subprogram  once  for  each  autoevent  in  the  network.  The  charac¬ 
teristics  of  these  subprograms  are  summarized  in  Table  6. 1. 

The  state  area  must  contain  three  items  •  the  state  set  (in  its 
multidimensional  form),  the  ring  whicii  identifies  the  state  variables, 
and  the  state  mapping  structure.  The  first  and  third  are  required 
for  the  state  mapper  (which  provides  a  mapped  state  when  a  vector  is 
given),  the  second  is  required  in  the  specificatien  of  results  (cf. ,  sec¬ 
tion  7).  All  three  are  required  in  the  generation  of  the  matrix  struc¬ 
ture  and  the  compilation  of  the  results  matrix.  Flgur®  6.  3  illustrates 
a  reasorable  form  for  the  structure  in  this  area.  The  symbol  table  is 
used  precisely  as  it  was  ir  the  retwerk  a.'^d  working  •  etv/ork  structures. 
The  procedure  for  constructing  *he  state  mapping  struc';ure  should  be 
clear  from  eqs.  6.  3  -  6.  6  . 

Two  difficulties  will  be  encountered  in  programming  ihe  matrix- 
event  generator.  In  converting  an  nver.t 

^  k  ^  ^k’  ‘^k 

to  matrix  form,  we  must  note  th«*  Trst.  although  a  part  of  bj^  may  be 
rectangular,  it  will  map  into  ’’block"  "j.  general  orJy  If  it  is  also  a 
subset  of  the  state  rectangles:  and,  second,  the  ’increment"  g^^  will 
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map  into  a  constant  only  when  it  does  not  represent  a  crossing  of  a 
boundary  between  state  rectangles.  Thus,  if  the  state  set  were  oi 
tv'o  dimensions  and  consisted  of  two  rectangles  Sj  and  Sg.  and  bj^ 
were  a  rectangle  contained  in  Sj.  bj^  may  still  need  to  be  represented 
by  more  than  a  single  ’’block"  in  the  matrix  as  illustrated  in  Figure 
6.  4.  Here  the  event  is  assumed  to  take  each  state  in  the  direction  of 
Sg- 


Sz 


Figure  6.  4 

Illustrating  "splitting"  of  bj^  for  constancy  of  g. 

All  but  the  rightmost  column  woulo  have  a  constant  mapped  incre¬ 
ment  gj^.  found  by  simply  weighting  the  i-th  component  of  gj^  by  the 
respective  6  j..  However,  it  can  be  readily  shown  that  the  mapped 
increment  (g'j^)  would  generally  be  different  for  each  of  the  states  in 
the  rightmost  column  of  b^^. 

Some  procedure  like  ttwt  diagrammed  in  Figure  6.  5  would  have 
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Figure  6.  5  Tu::.  matrix-eveni  generator. 
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to  be  used  to  take  account  of  these  difficulties.  Here,  one  first  deter¬ 
mines  which  of  the  rectangles  contains  which  portion  of  bj^.  These 
are  designated  b',  and  are  removed  from  the  b^^.  Then  one  determines 
which  of  the  b'  cross  into  a  different  state  rectangle.  These  are  desig¬ 
nated  b”,  and  consist  of  a  rectangle  or  set  of  rectangles  readily  mapped 
by  a  block  or  blocks.  The  remainder  of  the  b'  must  be  treated  indiv¬ 
idually  by  a  routine  (labeled  "incorporate  crossing  . . ,  "). 

Experience  with  handworked  examples  suggests  that  the  boundary 
crossing  cases  can  fairly  frequently  be  added  to  existing  blocks,  but 
in  a  nonpredictabie  way,  and  that  this  is  sufficiently  frequent  to  warrant 
considering  a  scan  in  the  "incorporate  crossing. . .  "  routine  which 
adds  stray  terms  to  existing  blocks.  However,  the  cost  of  such  a  scan 
and  the  typicalness  of  the  handworked  examples  are  still  imponderables. 
In  ary  event,  the  advantage  of  treating  transitions  within  a  state  rec¬ 
tangle  by  blocks  in  the  structure  "athcr  than  individually,  seems  un¬ 
questionable.  There  should  be  a  sm^ll  number  of  state  rectangles  rep¬ 
resenting  large  numbers  of  states,  giving  major  advantage  to  the  tech¬ 
nique. 

It  should  be  declared  also  that  since  execution  of  the  matrix- 
vector  multiplications  will  probably  rot  use  every  indexing  line  as  a 
register  (hence  fast)  operation,  and  since  order  of  indexing- lines  within 
a  block  is  irrelevant,  tl^  procedvn*es  should  pl^e_tjie  irdexing-Iine 
with  highest  repetitivity  'fi)  first  in  every  case,  on  the  assumption  that 
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only  the  first  one  will  involve  exclusively  register-to-register  indexing. 
The  extension  of  this  idea  if  more  registers  :\re  available  is  obvious. 

Upon  completion  of  the  execution  of  the  state  and  matrix  construc¬ 
tor,  the  working  nebvork  area  contains  no  useful  information,  and  will 
be  cleared. 
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7.  THE  SOLUTION  FOR  EQUILIBRIUM  PROBABILITIES 

Once  the  matrix  structure  has  been  formed,  it  can  be  used  to 
calculate  the  equilibrium  probability  of  each  state.  A  single  command 
will  initiate  the  process,  which  makes  use  of  thrtBe  areas :  the  matrix 
area,  a  state  probability  area,  and  a  working  probability  area.  The 
latter  two  areas,  which  will  contain,  vectors  of  state  probabilities, 
are  created  upon  entry  to  this  process.  The  workirig  probability  area 
is  temporary  and  will  be  destroyed  upon  completion  of  the  process. 

The  input  of  this  program  »s  the  matrix,  a  specification,  of  de¬ 
sired  accuracy,  and  a  specification  of  the  maximum  number  of  itera¬ 
tions  to  be  allowed  (which  provides  a  limitation  to  the  investment  in  a 
run,  protecting  the  user  from  excessive  expense  for  ill-concitioned 
problems\  The  output  is  the  equilibrium  state  probability  vector, 
which  is  left  in  the  state  probabri  ty  jrea  ard  an  error  estimate,  also 
placed  in  the  state  probability  Table  7.  1  summarizes  the  char¬ 

acteristics  of  this  program. 

The  numerical  technique  to  be  followed  ’S  an  iteration  process 
of  the  form 

:7  ,  :t  A  (7. 1 

nal  n  ' 

where  A  is  a  stochastic  matrix  ard  r 's  a  state  probability  vector. 
When  A  is  formed  from  a  matrix  of  trarsHion-intensities  in  a  partic¬ 
ular  way,  and  when  the  initial  Uerate  is  suitably  chosen,  then  the 
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limit  of  this  iteration  is  the  equilibrium  probability  vector 

•-  lim  tt  (7. 2) 

n  ' 

n-*oo 

A  detailed  treatment  of  the  numerical  technique  is  provided  in  [l], 

which  describes  the  prototype  program,  RQA-1. 

The  matrix  stored  in  the  matrix  area  is  not  precisely  the  matrix 

of  transition  intensities,  since  it  does  not  include  diagonal  elements. 

T 

Let  d  -  <dj, . . . ,  djj  >  be  a  column  vector  consisting  of  the  diag- 

S 

onal  elements.  Let  Q  be  the  matrix  actually  stored  in  *he  matrix 
area.  Then  the  matrix  U  of  transitional  intensities  is 

U  Q  -  D  (7.  3) 

where 

d  Q1  (7.4) 

and  1  is  a  column  vector  of  appropriate  dimension  whose  elements 
are  ail  unity.  Then,  according  to  eq.  [12'  of  [l],  with  a  suitable 
change  of  notation,  the  stochastic  matrix  A 's  given  by 

A  MQ  .  D}  i  1  (7.5) 

where 
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max  d. 


(7. 6) 


Thus  the  iteration  process  is 

7,-7  <  AQ'  -  7  AD  •*  -  . 
ns  1  r  ^  n  n 


(7.  7) 
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no 


This  program  follows  the  basic  outline  of  the  RQA-i  program, 
which  we  v/ill  assume  is  familiar  to  the  re^ider.  It  differs,  however, 
in  several  maicr  ways. 

First,  because  the  matrix  stored  describes  Q  rather  than  U, 
the  iteration  process  is  different.  The  values  of  A  and  d  must  be  cal¬ 
culated,  AQ  is  formed,  the  previous  iterate  :r^  is  copied  from  the 
probability  area  into  the  working  probability  area,  and  then  ff^(AQ) 
is  accumulated  onto  the  probability  area,  forming  :Tj.{AQ)  +  .  Fin¬ 

ally,  d  is  used  to  accumulate  onto  the  same  area,  to  complete 
eq.  (7.7). 

Second,  since  the  matrix  Q  is  formed  automatically,  there  is 
considerably  less  dema.'iu  for  testing  and  diagnostics,  and  the  sizes 
of  the  arrays  known  and  dc  not  reed  to  je  calculated. 

Third,  there  will  be  ro  repeated  runs  in  any  particular  sys¬ 
tematic  form,  so  thit  the  choice  of  how  tc  form  the  'r.itial  iterate 
is  narrowed.  It  is  also  r  arrowed  by  the  greater  remoteness  of  the 
user,  since  he  has  no  kr;Owledge  of  the  state  space  md  is  unlikely  to 
want  to  provide  his  own  initial  Iterate. 
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A  flow  diagram  for  this  program  is  shown  in  Figure  7. 1,  where  the 
names  of  corresponding  RQA-1  programs  are  inserted  (in  parentheses) 
where  appropriate.  A  default  is  provided  so  that,  if  a  matrix  has  not 
been  prepared,  the  program  supplies  a  trivial  solution.  This  prevents 
system  shutdown  over  such  errors,  with  the  error  having  a  clear  mean¬ 
ing  to  the  user. 

It  is  expected  that  an  interface  program  will  be  created  which 
will  allow  a  batch  or  teletype  user  to  directly  supply  a  list  of  auto¬ 
events  and  state  rectangles,  and  a  description  of  a  results  matrix  (see 
section  8),  and  obtain  a  solution  without  involving  the  compiler  and  the 
network  generation  programs.  (In  this  way  we  obtain  a  very  powerful, 
primitive  model  solver  which  is  akin  to,  but  an  improved  version  of, 
RQA-1. ;  This  interface  program,  which  calls  on  several  QAS  sub¬ 
routines  including  the  "State  Probability  Solver, "  will  be  called  RQA-2. 
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8.  SPECIFICATION  AND  CALCULATION  OF  RESULTS 


The  results  desired  by  the  user  of  this  system  will  not  ordinarily 
be  the  state  probabilities  ir  their  raw  form.  Rather,  they  will  be  vari¬ 
ous  quantities  which  are  readily  computed  from  the  state  probabilities. 
For  example,  one  may  desire  the  probabilHy  distribution  of  the  number 
of  tasks  in  a  particular  queue;  or  the  mean  queue  length;  or  the  proba¬ 
bility  that  a  particular  server  is  ir  use,  or  the  mean  rate  of  service 
completions  in  that  server,  or  the  rate  of  transfer  of  tasks  through  a 
particular  port;  or  the  probabiUty  that  a  task  is  blocked  at  that  port. 

Each  of  these  represents  a  weighted  sum  of  state  probabilities, 
or  a  set  of  weighted  sums.  Tr  other  words,  there  is  a  matrix  such 
that  the  result  vector  w  'perhaps  dimension  one)  is 

w  -  r*"  .  (8.1) 

!r  addition  many  other  useful  »^esu}*s  whU^  ro^  of  this  tvpe  are  »'e<td- 
Uy  found  from  results  W’hich  *re  of  type.  For  example,  the  ex¬ 
pected  wait’pg  t'me  m  a  queue  wbVh  is  the  mean  rumbt  r  of  tasks  ir. 
the  queue,  divided  by  the  mean  rate  of  t,sk  transfer  Tto  its  ’nput); 
or  the  mear  time  between  service  ccmrle^iors  of  a  server  '.which  is 
the  reciprocal  of  the  mean  rate  of  service  rcmpletiors).  Such  results 
can  be  treated  as  .functions  of  weighted  sums 

W  -  fiw)  (8.2) 

which  are  desc’*ibable  by  subrcut’nes. 
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Fi^rthermore,  in  any  particular  solution,  the  user  may  have  sev¬ 
eral  results  which  he  would  like  to  observe  simultaneously,  so  that  the 
vector  w  may  consist  of  several  results  strung  together,  and  the  f;;nc- 
tion  f  may  be  a  set  of  functions  on  different  parts  of  the  vector. 

There  are  four  separate  operations  which  QAS  must  be  prepared 
to  carry  out  to  produce  output: 

(1)  The  specifications  of  results  must  be  retained  in  a  storage 
structure. 

(2)  The  matrix  F  must  be  pre{»red. 

(3)  The  vector  w  must  be  computed. 

(4)  The  result  w  must  be  transformed  into  the  format  required 
for  SELMA. 

The  commands  associated  with  the  first  operation  will  be  added  to  the 
generation  phase  commands,  while  the  last  three  of  these  operations 
will  together  constitute  a  phase  called  the  results  rfiase. 

^  The  Spec^f^cat^or  of  Results 

The  roed  for  a  separate  structure  which  "remembers”  the  re¬ 
quested  results  stems  primarily  from  the  desirabilHy  of  allowing 
SELMA  users  to  specify  results  at  any  time  during  the  development 
of  the  network  model.  Often  he  will  want  comparable  results  for  sev¬ 
eral  modifications  of  the  network,  and  will  want  to  change  the  network 
without  charg’r.g  specif icatiors.  Often,  too,  he  will  want  to  make  minor 
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alterations,  as  an  afterthought,  long  after  he  has  decided  v^hat  results 
are  to  be  observed.  The  specification  can  not  be  transformed  to  its 
final  form  (a  matrix  '  )  until  the  network  is  complete  and  the  state 

mapping  structure  determined.  Thus,  the  results  specifications  must 
be  saved,  awaiting  a  command  to  prepare  the  matrix  and  to  solve  for 
the  results  at  the  appropriate  time. 

The  structure  which  accomplishes  this  shou  ivailable  during 
the  generation  of  the  network,  however,  and  const  _ntiy  should  be 
placed  in  the  network  area.  During  generation,  changes  in  the  net¬ 
work  may  invalidate  existing  result  specifications,  and  procedures 
which  at  e  responsible  should  make  adjustments.  For  example,  a  queue 
whose  expected  length  has  been  specified  as  a  result  may  subsequently 
be  eliminated  from  the  network.  During  the  operation  of  the  ’’destroy 
element”  command  in  network  generation  phase,  a  test  should  be  maae 
to  determine  if  a  result  specification  is  involved,  and  if  so  to  eliminate 
the  result  specification  at  the  same  time. 

T-he  reference  of  a  specification  to  an  element  or  a  port  is  syn¬ 
tactically  like  a  connection,  in  that  a  aestruction  of  tne  object  requires 
a  choice  of  destroying  the  connection  or  refusing  tc  oestroy  the  object. 
Thus,  one  form  which  this  structure  can  take  is  a  ring  structure  form, 
aading  new  attributes  to  the  element  ana‘'or  port  blocks  ana  connecting 
them  to  appropriate  ’’result  blocks.  ”  Since  the  number  of  results 
specifieu  will  generally  be  very  .small,  most  such  attributes  would 
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normally  have  their  default  (empty)  values.  A  more  efficient,  but 
less  uniform,  structure  would  provide  a  simple  table  of  result  specifi¬ 
cations  which  can  be  searched  by  generation  routines  ana  rewritten 
in  the  (improbable)  event  that  alterations  are  required. 

A  minimum  set  of  result  specifications  is  tabulated  in  Table  8. 1. 
The  purpose  of  each  of  these  specifications  shoula  be  self-evident, 
except  for  the  sixth.  By  means  of  this  latter  specification,  it  is  pos¬ 
sible  to  establish  the  probability  that  an  input  port  is  in  a  blocking  state, 
or  that  an  output  port  has  a  blocked  task.  All  of  the  specifications  in 
the  table  represent  weighted  sums  of  state  probabilities. 

Since  commands  altering  the  results  specifications  will  be  inter¬ 
mixed  arbitrarily  with  generation  commands,  and  their  structures  are 
in  the  same  area,  these  commands  are  best  treated  as  additional  rou¬ 
tines  in  the  generation  phase  program.  Table  8.  2  inaicates  a  minimal 
set  of  results  specification  routines  needed  to  carry  out  alterations 
of  the  results  specifications. 

8.  2  The  Calculation  of  Results 

The  calculation  of  results,  which  can  result  from  a  single  QAS 
command  issued  after  tne  compilation  of  a  network,  involves,  as  we've 
said,  three  operations:  preparation  of  the  matrix  .  computation 
of  the  vector  w  (cf.,  eq.  8. 1  ).  and  formatting  for  output  to  SELMA. 

In  this  section,  we  will  assume  that  the  last  operation  constitutes  a 
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Table  8. 1  Result  Specifications 
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simple  supplying  of  the  appropriate  portions  of  the  vector  w  with  the 
copies  of  the  results  specifications  so  that  the  result  can  be  identified 
by  SELMA.  Thus>  we  assume  that  SELMA  takes  responsibility  for  mak¬ 
ing  necessary  further  transformations,  either  of  the  form  of  eq.  (8.  2) 
or  into  graphical  displayable  form. 

(If  this  is  not  the  case,  and  if  QAS  is  to  provide  more  general  re¬ 
sults  than  those  which  were  described  in  Table  8.  2,  or  ii  QAS  must  pre¬ 
pare  the  output  graphics,  then  further  structures  must  be  supplied  to 
accept  format  specifications  and  function  identifications  in  the  results 
calculation  phase. ) 

The  commands  for  this  phase  are  summarized  in  Table  8.  3.  The 
only  nontrivial  operation  in  the  figure  is  the  preparation  of  the 
matrix  T  .  For  this  purpose  an  area  is  created,  called  the  results 
matrix  area,  where  the  structure  describirg  !  will  be  kept.  The 
matrix  f  should  probably  be  represented  ir.  matrix  outline  form,  with 
the  modification  that  the  increments  g  should  be  replaced  simply  by  the 
row  index  of  the  result.  This  reflects  the  fact  that  repetitive  element 
values  in  T  will  generally  be  along  columrs  instead  of  diagonals.  Nat¬ 
urally,  the  vector-matrix  multiply  operator  will  reed  tc  be  different , 
reflecting  this  change. 

To  illustrate  how  this  process  would  proceed,  consider  a  column 
of  the  matrix  corresponding  to  the  probibiTty  that  a  particular  queue 
has  length  two.  The  elements  of  ♦h's  column  will  consist  of  ones  at 
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each  row  corresponding  lO  a  state  for  which  the  queue  length  is  two. 

This  set  of  row  indices  is  the  set  found  by  restricting  the  state  set 
the  value  2  in  the  appropriate  i»niension.  An  operator  ("set  restric¬ 
tion")  for  creating  such  a  set  of  n-tuples  already  exists  (see  Table 
4. 4),  and  is  readily  used.  A  set  of  iteration  blocks  for  the  matrix 
outline  corresponding  to  this  column  can  be  generated  from  the  set 
of  n-tuples  using  programs  developed  foi  ihe  transition  matrix  trans¬ 
lation. 

The  operation  of  preparing  columns  of  the  results -matrix  are 
correspondingly  simple  for  every  specification-type  in  Table  8. 1  ex¬ 
cept  type  5.  For  this  type,  the  flow  rate  through  a  port  is  to  be  cal¬ 
culated.  Equivalently,  the  probability  intensity  of  occurrence  of  a 
particular  spontaneous  event  must  be  shown  for  each  state  of  the  net¬ 
work.  Since  the  autoevents  which  can  cause  the  spontaneous  event 
are  not  obvious,  the  desired  probability  intensities  are  not  easily 
found,  particularly  from  the  information  in  the  network  structure. 

Nevertheless,  the  result  specified  in  type  5  specification  is  val¬ 
uable,  and  worthy  of  the  difficulties.  What  is  required  is  a  form  of  par¬ 
tial  compilation  of  the  network,  exploring  only  those  connections  whose 
spontaneous  events  can  cause  the  spontaneous  event  under  examination. 
By  this  procedure,  each  autoevent  of  the  associates  of  the  connection 
are  examined  to  find  the  subset  of  their  autocondition  sets  which  are 
"taken  into"  the  forbidden  states  of  the  connection  by  the  autoevent. 
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The  intensity  (fi)  of  this  autoevent  is  aided  to  the  matrix  outline  for  the 
state  subset  so  found.  The  exoevents  of  tne  associates  are  similarly 
examined,  and  whenever  a  nonempty  subset  of  the  endocondition  set  at 
a  port  is  found  which  is  ’’taken  into"  the  forbiaden  sei  by  the  exoevent, 
the  connection  at  the  offending  port  is  absorbed.  The  new  autoevents 
are  then  examined  as  before,  and  the  process  repeals  until  no  more 
absorptions  are  necessary.  This  process  will  usually  involve  absorp¬ 
tion  of  only  one  or  two  of  the  connections.  £  nd  can  make  g(iod  use  of 
existing  compiler  subprograms.  All  the  intermediate  results  can  be 
kept  in  the  working  network  area  *s  the  compiler  does. 

Alternatively,  the  network  can  be  compiled  with  the  connection  in 
question  being  absorbed  last.  The  desired  information  is  readily 
found  in  the  process  of  the  final  absorption,  for  which  a  special  sub¬ 
routine  can  be  written  that  prepares  the  column  of  ■  instead  of  act‘~;aiiy 
absorbing  the  connectioa  Of  course,  this  method  is  expensive  since  it 
duplicates  the  process  of  compilation. 

The  program  which  prepares  the  results  matrix  continues  to  add¬ 
on  columns  to  the  f  matrix,  suitably  identified  with  the  specifications, 
until  all  specifications  have  been  treated.  Thus  all  results  are  computed 
simultaneously  when  the  matrix  is  multiplied  by  the  equilibrium  probabil¬ 
ity  vector. 

The  computation  of  the  matrix  as  an  intermediate  result,  rather 
than  computing  the  results  directly,  will  facilitate  future  generalization 
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of  QAS  whereby  parameter  values  will  be  permitted  to  be  changed  without 
necessitating  recompilation  of  the  transition-matrix  or  the  results- 
matrix. 


9.  EPI1.OGUE 

It  has  been  necessary  to  describe  the  system  in  considerable  detail 
in  order  to  work  out  the  critical  interrelationships  between  the  major  pro¬ 
cesses  and  data  structures.  In  this  concluding  section,  a  series  of  ob¬ 
servations  will  be  made  about  some  auxiliary  matters  which  are  not  cen¬ 
tral  to  the  main  interrelationships,  and  hence  do  not  need  to  he  presented 
in  such  detail.  These  observations  represent  a  summary  of  organization 
concepts,  notions  for  useful  generalization,  and  certain  reservations  about 
unresolved  theoretical  questions. 

9. 1  Organization 

We  have  described  a  system  for  the  solution  of  simplo  stochastic  net¬ 
works  as  an  amorphous  collection  of  well-defined  routines  which  operate 
upon  collection  of  well  defined  data  structures.  The  organization  of 
this  system  is  expected  to  be  imposed  by  the  calling  system  (in  this  case 
SELMA),  and  ultimately  subject  to  the  whim  of  a  user.  Such  an  organiza¬ 
tional  philosoi^y  is  a  natural  one  for  a  progivimming  system  designed 
for  conversational  use,  and  imposes  a  minimum  constraint  upon  the 
thought  sequence  of  the  user. 

To  implement  this  organizational  philosophy,  each  routine  should 
be  designed  to  function  in  any  position  in  a  command  sequence  without 
catastrophe.  The  consequence  of  the  routine's  operation  should  be  the 
most  logical  and  consistent  one  that  could  be  expected.  In  some  cases 
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that  m’ght  htt’'?  tc  be  ar  error  diagnostic,  ard  x  def^ul*  operation.  In 
ar.y  case,  the  "esp^'cse  must  be  corsister.t  w'*h  pre’'’cus  operarcns. 

For  e.vampJe  'f  i  ’  gereratior  commix'd''  [altering  the  network) 
is  »'ece  ved  afie'*  x  <"-'mr'le  ccmm.xnd.  ard  fher.  ri.  "CAlciIate-prob- 
abilities’  ccmm.*rd  ^  received  the  most  logical  response  to  the  latter 
is  to  or'^v’de  the  probabilities  for  the  altered  network,  rather  than  use 
the  ’'esuUs  of  *he  ‘V''mr  le  ‘  ''t**ar.s*t»C’'  mitrix"'.  Thus,  upon  alter¬ 

ing  a  retwc-k  trarsiticr  matrix  fr-^m  ary  previous  ’’compile” 
should  be  des’rcyed  ard  the  'c alcuUte-prc  b:4b:l:t:es”  should  either 
g:v-’  a'-  errrr  u  d«r..*t’o'^  [ror  fatal[  when,  ro  trarsH''^r  matrix  is  pre- 
sen*  O’-  sh'  u'd  a  i*‘'>m<i*.'callv  call  the  ”comp:’e”  program.  (If,  for 
•'‘>aS''r  th  :  'jser  des*''es  tho,*  ♦he  t’':r.sH;or  mat-ix  be  saved  be- 
fc’'e  the  re**vo’'k ’s  tl*e’’ed  h'*  should  use  a  ’'drc  imertatlor”  command. ) 

S  m  enter  *•■/=,}  rre ’■s's‘err-'.es  error  ir  t^e  alte^aticr  of  results 
spec*'*rj=''  q  c -Irul  .t*''"  of  resu.’*s  d  the  use  of  sU^-te  mapping 
st»'uc tune  ^r.i  pnob^b’'/*ies.  All  '‘f  ♦Jicse  rrrflic*s  c-r  be  re¬ 
solved  by  ’rs*s’  'g  th.t  ’he  nrfcrraa*’or  t  a!!  x"-'ix5  be  appVc.ble  at 
ail  t'm^s.  tn  the  c  inrer*  network  *rd  results  srr’C’f'c  ard  auto- 
m.:.t  c  j.ltv  rln_r>  g  such  4,’'eas  v.-h^»- jver  ’rcors'ste’ry  ^r’ses.  !f 
this 's  dor  e,  no  ccrtrol  over  the  order  of  comm.xrds  will  be  ’’equired 
of  SELM.-1  by  Q/^S.  ,SELMA  miV  howe^e-  mpose  such  rortrol  H- 
self.  [ 

Fr'r  (  c' ver’-^rre  the  romm.i’ds  and  jssoc’lated  irfiases  and 
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areas  are  I’Sted  in  Table.  9.  1.  In  the  column  marked  ’’areas  cleared,  ” 
areas  in  parentheses  are  those  automatically  cleared  because  of  in¬ 
duced  inconsistency. 

9.  2  Documentation 

Documentation  phase  commands  will  provide  the  ability  of  a  user 
to  save  any  network  structure,  transition  matrix-state  mapping,  results 
specification,  or  results  as  a  named  file,  and  to  restore  them  to  QAS 
for  resumption  of  analysis.  This  documentation  must  be  done  in  such 
a  way  that  QAS  car,  at  a  later  time,  be  restored  to  an  operating  condi¬ 
tion  with  the  saved  network. 

Moreover,  the  file  will  have  to  contain  sufficient  iriormation  to 
restore  the  displays  of  SELMA,  via  dataphone.  since  the  PDP-9  has  no 
facility  for  long-term  saving  of  display  files. 

9.  3  Deferred  Evaluation  of  Par;am^ters 

One  of  the  most  cruc’al  imprcvemet  is  which  this  system  require.e 
is  a  provision  which  permits  compilaticr  of  networks  and  results  speci¬ 
fication  with  parameters  in  the  elements'  parameter  lists  treated  as 
variables.  With  such  a  prcvis’on,  the  operation  rf  compilation  will 
be  expanded  'rto  thi’ee  operations:  compilation,  evaluation,  and  calculation. 
The  first  operation,  compiles  the  network,  preserving  an  arithmetic  des¬ 
cription  of  how  each  parameter -dooendert  constant  can  be  found.  The 
second  operafior  uses  the  arithmetic  descriptions  to  place  the  parameter- 
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dependent  constants  in  the  appropriate  locations  of  the  transition  matrix 
structure.  Then,  when  a  user  is  exploring  a  single  model  with  many 
values  of  the  parameters,  the  compilation  is  performed  only  once  and 
the  evaluation  (and  calculations)  are  performed  for  each  set  of  parameter 
values. 

With  anticipated  modes  of  use,  based  upon  current  usage  of  RQA-1 
which  has  a  similar  feature,  this  will  reduce  the  number  of  compilations 
required  by  an  order  of  magnitude  or  more.  Since  compilation  is  an  ex¬ 
tensive  operation,  this  is  an  important  economic  improvement,  and  will 
also  provide  a  dramatic  improvement  of  response  time. 

Indications  are  that  compile  time  will  not  be  dramatically  in¬ 
creased  by  implementing  this  capability,  and  that  evaluation  can  take 
considerably  less  time  than  compilation.  An  experimental  prototype 
of  one  system  for  accomplishing  this  was  developed  [8]  for  use  on  the 
IBM  7090  and  demonstrated  feasibility. 


REFERENCES 


1.  Wallace,  V.  L.  and  R.  S.  Rosenberg,  RQA-1,  The  Reciu^sive 
Queue  Analyzer,  Technical  Report  2,  Systems  Engineeriug  Lab¬ 
oratory,  Department  of  Electrical  Engineering,  University  of 
Michigan,  Ann  Arbor,  February  1966. 

2.  Wallace,  V.  L.  and  R.  S.  Rosenberg,  "Markovian  Models  and 
Numerical  Analysis  of  Computer  System  Behavior,"  proc. 

IFIPS  Spring  Joint  Computer  Conference,  vol.  28,  1966,  pp. 
141-148. 

3.  .Jackson,  J.  H. ,  "Matrix  Storage  Scheme  Suitable  for  SysUim 
360, "  Memo  to  07449  File,  Systems  Engineering  Laboratory, 
University  of  Michigan,  Ann  Arbor,  December  3,  1965. 

4.  MTS,  Michigan  Terminal  System,  Computing  Center,  Univer¬ 
sity  of  Michigan,  Second  Edition,  December  1,  1967. 

5.  Newman,  W.  M. ,  The  ASP-7  Ring  Structure  Processor,  Com¬ 
puter  Technology  Group  Report  67/8.  Imperial  College,  London, 
October  1967. 

6.  Wallace,  V.  L. ,  "Preliminar7  Description  of  a  Free  Core  Mana¬ 
ger,  "  Memo  to  07842  File,  Systems  Engineermg  Laboratory, 
University  of  Michigan,  Ann  Arbor,  November  9,  1967. 

7.  Sutherland,  W.  R. ,  The  On-Line  Graphical  Spet  ifications  of  Ccr  - 
puter  Procedures,  Lincoln  Laboratory  Technical  Report  405,.  " 
M,  I.  T- ,  Lexington,  Massachusetts ,  May  1966.  (iippendix  C, 
"The  CORAL  Language  and  Data  Structure"). 

8.  Jackson,  J.  H. ,  "Internal  Structure  of  the  IBM  7090  Experimental 
Version  of  the  Deferred  Evaluation  System  (DES-0), "  Memo  to 
07842  File,  Systems  Engineering  Laboratory,  The  University  of 
Michigan,  Ann  Arbor,  September  15,  1966. 

9.  Randall,  L.  S. ,  I.  S.  Uppal,  G.  A.  McClain,  J.  F.  Blinn,  An 
Implementation  of  the  Queue  Analyzer  System  (QAS)  on  the  IBM 
360/67,  Systems  Engineering  Laboratory  Technical  Report  43, 
Department  of  Electrical  Engineer ing,  University  of  Michigan, 

Ann  Arbor,  1969.  Concomp  Project  Technical  Report  22.  Com¬ 
puting  Center,  University  of  Michigan,  Ann  Arbor,  1969. 


131 


132 


iO,  WaUice,  V.  L. ,  On  the  Represer-tattcr.  of  Markcv’ir.  Systems  by 
Network  Models  Systems  Engineering  Laboratory  Tec hrJcai 
Reiiort  ¥i~,  Department  of  Electrical  Engineering,  Un’versity  of 
Michigan,  Ann  ArboT*.  1939;  also  Ccnccmp  Project  Technical 
Report  21.  Ccmputlrg  Center,  UnH'ersi^v  of  Michigan,  An*- Arbor, 
1939. 


11.  Jackson.,  J.  H. ,  SELMA:  A  Conversat'cnrJ  System  for  Graphic  A» 
Specification  of  M^kcviar.  Queueing  Networks ,  Systems”  EngineeH'' g 
Labor  Jory  Technicitl  Report  45,  Department  of  Electrical  Er.gi- 
reerirg.  University  of  Michigan,  Ann  Arl)or,  1969,  also  Concomp 
Technical  Report  ?3.  Comput’ng  Center,  Urh'e~S’ty  ofM;r.h-gan, 

Ann  Arbor,  1939. 


APPENDIX  A  . 


Preliminary  Description  of  a  Free-Core  Manager 

This  is  a  brief  proposal  for  a  programming  system  which  gives 
and  collects  blocks  of  free  core  (hiring  execution  of  list- manipulating 
programs.  It  is  assumed  here  that  three  properties  are  required  by  the 
programs  which  use  this  system: 

1.  Blocks  of  core  in  use  must  be  format  free.  In  other 
words,  there  may  be  no  reserv*  d  words  or  bits  in 
blocks  which  are  in  use. 

2.  Blocks  of  arbitrary  size  may  be  requested  or  released, 
but  certain  sizes  will  be  much  more  frequently  used  than 
others. 

3.  No  relocation  is  possible  for  a  block  once  it  is  put  to  use. 

This  proposed  system  is  believed  to  be  a  satisfactory  compromise 

between  time,  storage  efficiency,  and  programming  effort  under  the 
above  circumstances. 

The  area  in  which  the  blocks  are  to  be  allocated  will  be  con¬ 
sidered  a  large,  named  vector.  Thus,  if  there  is  reason  to  segregate 
blocks  for  economic  or  other  reasons  to  different  parts  of  virtual 
memory,  it  is  possible  to  exert  such  control.  (For  example:  different 
parts  of  the  list  structures  to  be  generated  may  be  in  use  at  different 
times,  and  it  may  be  uneconomical  to  keep  them  all  in  virtual  memory 
at  once. ) 
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The  first  n  words  of  each  area  should  represent  some  number  of 
categories  (say  eight) .  Half  of  each  of  these  words  will  be  pointers  to 
a  member  of  the  category  or  back  to  itself  if  there  are  no  members. 
That  is,  they  each  form  the  head  of  a  ring.  The  other  half  word 
indicates  th  number  of  members  in  the  rii^.  The  members  of  the 
categories  are  blocks  of  free  memory  in  a  particular  range  of  sizes. 
(Usually,  most  categories  will  consist  of  a  single  commonly  used 
size,  with  one  category  being  a  '’miscellaneous"  category. )  These 
free  memory  blocks  will  have  the  format  shown  in  Figure  A  1.  Notice 
that  one-word  "blocks"  have  a  special  format.  For  this  reason,  it  is 
necessary  that  one-word  blocks  always  be  in  the  first  category. 

There  should  be  at  least  three  calls  to  the  manager  system 

ALLOC  (AREA,  SIZE) 

GET  (SIZE,  AREA.  ERROR) 

FREE  (LOC,  SIZE,  ERROR) 

Their  functions  will  be  as  follows: 

ALLOC  (AREA,  SIZE) 

Here  "AREA"  is  the  name  of  the  vector  of  length  SIZE  which  is 
to  be  used  by  the  manager  as  its  working  region,  SIZE  must  be  less 
than  or  equal  to  the  dimension  declared  (elsewhere,  in  the  using 
program)  for  the  vector.  ALLOC  may  be  used  as  often  as  desired, 
setting  up  as  many  regions  as  necessary.  ALLOC  sets  up  the  category 
heads  in  the  vector  suitably  linked  so  that  all  categories  are  empty 
except  "misc"  which  contains  one  large  free  block  representing  the 
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rest  of  the  vector.  Of  course,  ALLOC  also  puts  the  right  head  and  tail 
words  on  this  free  block. 

GET  (SIZE,  AREA,  ERROR) 

GET  is  a  function  which  locates  a  free  block  of  the  specified  size 
in  the  specified  area.  It  has  value  /-qual  to  the  address  of  the  first 
word  in  the  located  block.  If  a  block  of  the  exact  size  is  not  available, 
GET  splits  a  larger  block.  If  no  larger  block  is  available  a  return  is 
made  to  the  specified  error  statement  label.  The  ond  and  third 
arguments  are  optional  with  default  values  "same  as  last  specified 
in  an  ALLOC,  GET  or  FREE  call",  and  "system  return",  respectively. 

The  procedure  for  GET  is  charted  in  Figure  A  2.  The  operations 
marked  with  an  asterisk  (*)  are  adjustable  procedures  which  can  be 
controlled  so  that  the  sizes  of  remainders  after  splits  have  a  dis> 
tribution  close  to  that  of  expected  requests.  They  should  be  under 
some  control  of  the  user,  as  also  the  choice  of  sizes  in  each  category 
(which  is  made  by  a  procedure  in  the  FREE  subroutine) . 

FREE  (LOC,  SIZE,  ERROR) 

FREE  is  a  subroutine  which  places  a  block  back  into  an  approp¬ 
riate  category  of  free  core.  LOC  is  the  address  of  the  first  word  of 
the  block  to  be  released,  SIZE  is  its  size  (an  imeger) ,  and  ERROR 
is  the  statement  label  to  which  return  is  made  when  the  entire 
block  is  not  contained  in  a  current  working  region.  The  latter 
argumenc  is  optional,  with  default  value  "system  return". 
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In  order  to  avoid  splintering  of  free  core,  which  would  result  after 
a  long  sequence  of  GET's  and  FREE’S,  FREE  joins  the  block  being 
released  to  its  neighbors  if  they  are  also  free.  The  chart  of  Figure  A3 
shows  how  this  is  accomplished.  The  symbols  "Fl(. . . )  "  and  "F2(. . . ) " 
represent  the  left  half  word  and  the  right  half  word  of  ... ,  respectively. 

The  operation  marked  with  an  asterisk  {*)  has  a  portion,  loosely 
under  user  control,  which  decides  which  size  blocks  to  put  in  which 
categories.  It  is  through  the  "**’  programs  that  the  assignment 
strategies  can  be  adjusted  to  minimize  unnecessary  splintering  by 
taking  account  of  statistical  conditions  of  the  user's  list  structure. 

For  example,  use  of  categories  should  be  with  approximately  equal 
probabilities,  and  splits  of  large  blocks  to  assign  smaller  ones  should 
generally  produce  remainder  blocks  of  a  size  which  is  likely  to  be 
needed. 
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